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SQUARE FUNCTION ESTIMATES ON LAYER POTENTIALS 
FOR HIGHER-ORDER ELLIPTIC EQUATIONS 

ARIEL BARTON, STEVE HOFMANN, AND SVITLANA MAYBORODA 


Abstract. In this paper we establish square-function estimates on the double 
and single layer potentials for divergence-form elliptic operators, of arbitrary 
even order 2m, with variable t-independent coefficients in the upper half-space. 
This generalizes known results for variable-coefficient second-order operators, 
and also for constant-coefficient higher-order operators. 
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1. Introduction 

With this paper we continue towards the grand goal of resolving the Dirichlet 
and Neumann problems for general divergence-form higher order elliptic operators 
with LP data. The investigation of the second-order case has essentially spanned 
the past three decades in the subject, drawing from and giving back many cutting- 
edge tools of harmonic analysis, and by now the real coefficient case is relatively 
well understood. However, even the simplest higher order operators, such as the 
bilaplacian, presented outstanding difficulties: for instance, the sharp range of p 
for which the Dirichlet problem for the bilaplacian is well-posed in is still not 
known in high dimensions. And, to the best of our knowledge, no well-posedness 
results are currently available in the general variable coefficient case. 

In this work we aim to develop the method of layer potentials for general diver¬ 
gence form higher order elliptic operators. The main results of the present paper 
are square function estimates for single and double layer potentials in and the 
corresponding Sobolev spaces. We remark that one of the key difficulties in such a 
general context lies in the definition of suitable layer potentials and, more gener¬ 
ally, of Dirichlet and Neumann boundary data, as in the higher order case there is 
considerable ambiguity, some choices leading to ill-posed problems. Our approach 
is new even in the constant coefficient context, but is carefully crafted to handle 
the general case. 

Let us discuss the background and the results in more detail. 

In this project we study elliptic differential operators of the form 

( 1 . 1 ) Lu = {-ir Y. 

|a| = |/3|=m 

for m > 2, with general bounded measurable coefficients. As mentioned above, con¬ 
trary to the second order case, most of the known well-posedness results for higher 
order boundary value problems have been established only in the case of constant co¬ 
efficients (see, for example, |Ver90[ lPV95a[ rVer961 IShe06a[ ISheOBbl IKSlli IMM13b] . 
or the survey paper |BM15) b or concern boundary-value problems such as the 
Dirichlet problem 

(1.2) Lu=(-1)™ Y 5“(A„^5'^M) = 0inD, = / on SD 

\a\ = \p\=m 
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where the boundary data / lies in a fractional smoothness space. See |Agr07| 
IMMSlOl IBMlSb] . We are interested in the Dirichlet problem (|1.2I1 , with variable 
coefficients, in the classical case where the boundary data f lies in Li^idVl). 

1.1. The method of layer potentials, general framework. Classic tools for 
solving second-order boundary-value problems are the double and single layer po¬ 
tentials, given by 

(1.3) V^fiX)= [ vA*{Y)XYE^*{Y,X)f{Y)da{Y), 

Jon 

(1.4) S^g{X)= [ E\X,Y)g{Y)da{Y) 

Jan 

where u is the unit outward normal to 17 and where E^{X,Y) is the fundamental 
solution for the operator L. Making sense of formulas m and in our general 
context is one of the key tasks in its own right, and we will return to it below. 

It may be shown that, for any nice functions / and g defined on 517, the functions 
u = 'DqJ or It = S^g satisfy Lu = 0 away from 517. The classic method of layer 
potentials for solving boundary-value problems is to show that Sq (or Vq) satisfies 
some estimate, for example, 

[ |V5„+i5jf5'(A)pdist(A, 517) dX < C\\g\\L2(dQ), 

Jn 

and to show that the operator g i—>■ SQg\g^ is invertible from some space to another, 
for example, L'^{dfl) >->■ (517). Then the function u = ((5^ 1^^.^)“^/) is a 

solution to the Dirichlet regularity problem 

Lu = 0 in 17, u = / on 517, J |V5„+iu(X)p dist(X, 517) dX < (711/11^2(3^2). 

This method has been used in |Ver841 IDK871IFMM981 IZanOO] in the case of har¬ 
monic functions (that is, the case A = I and L = — A). This m ethod has also been 
used to study more general second order problems in [AAA+lfl IBarl3[ IBMl IRosl31 
IHKMPl3al IHMMlSl IHMM13] under various assumptions on the coefficients A. 
Layer potentials have been used in other ways in |KR091 IRul07[ IMit08[ |Agr09| 
IMMlll fAM14| . In particular, the second-order double and single layer potentials 
have been used to study higher-order differential equations in [PV921 IBMI3a] . 

1.2. Outline of the main results. In this paper we begin to generalize this 
method to the case of higher-order equations by defining the double and single 
layer potentials and for higher-order equations (see Section 1^ . and then 
by establishing some bounds on these potentials under certain conditions on the 
coefficients A. We plan to establish invertibility of layer potentials for some A, and 
thus well-posedness of the corresponding boundary value problems, in a forthcoming 
paper. 

Even in the case of second-order equations, some regularity assumption must 
be imposed on the coefficients A in order to ensure well-posedness of boundary- 
value problems. See the classic example of Caffarelli, Fabes, and Kenig [CFK81] . in 
which real, symmetric, bounded, continuous, elliptic coefficients A are constructed 
for which the Dirichlet problem with L^ boundary data is not well-posed in the 
unit disk. A common starting regularity condition is t-independence, that is, 

(1.5) A{x, t) = A{x, s) = A{x) for all x S K" and all s, t S R. 
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Boundary value problems for such coefficients have been investigated extensively in 
domains where the distinguished t-direction is always transverse to the boundary, 
that is, n = I(cc, t) : t > ip(a:)| for some Lipschitz function ip. See, for example, 

[JKMI IFJK841IKP931 IAAA+ 1 ll lAAHOSMAAMlOl \KKn[ lAR 121IAM141 IHKMPlbbl 

IHKMP l3al IBM] . (In two dimensions some well-posedness results are available even 
if the distinguished direction is not transverse to the boundary; see [KKPTpni 
IRuin7llB^?l3] .l 

The main result of this paper is the following theorem. 

Theorem 1.6. Suppose that L is an elliptic operator associated with coefficients 
A that are t-independent in the sense of formula (HI) and satisfy the ellipticity 
eonditions (j2.4ll and (1^ . 

Then the single and double layer potentials 5^ and T>^ in the half-space, as 
defined by formulas (12.321) and (12.351) . extend to operators that satisfy the bounds 

(1.7) f f \V^dtS^g{x,t)\^\t\dtdx<C\\g\\l2(^ny 

Jr^ J-cyo 

(1.8) [ r\V^dtV^n^,t)\^\t\dtdx<C\\f\\l^, (r)=qi/lli2(R„) 

Jr"J-oo ’ 

for all g € L^(R") and all f € |(®"); where C depends only on the dimension 

n + 1 and the ellipticity constants A and A in the bounds (i2a and (1^ . 

The space |(®") C L^(IR.") will be defined in Section 1^31 
We conjecture that this theorem may be generalized from the half-space to Lip¬ 
schitz graph domains, but the method of proof at the moment requires the extra 
structure of In the case of second-order operators (the case m = 1), bounds 

in the upper half-space may be immediately extended to bounds in domains above 
Lipschitz graphs via a change of variables, and so extra arguments are not neces¬ 
sary. In the higher-order case, this is not true, as the divergence form (HD is not 
preserved under changes of variables. (A different form of higher-order operator is 
preserved under changes of variables; such operators were investigated in |BM13a) .l 

1.3. Boundedness of layer potentials for second order elliptic operators. 

We now turn to the history of this problem. A reader familiar with the second 
order case may skip this subsection. As discussed above, layer potentials have been 
used extensively in the theory of second-order and constant coefficient boundary 
value problems. Thus, boundedness results for layer potentials have long been 
of interest. The celebrated result of Coifmann, McIntosh and Meyer [CMM82] 
established boundedness of the Cauchy integral on a Lipschitz curve; this implies 
that the operators / >->■ and 5 >->■ v • are bounded >->■ 

L^{dfl), where 11 is a Lipschitz domain and where A = J is the identity matrix 
(that is, where L = — A is the Laplace operator). From there many other bounds 
on harmonic layer potentials may be derived. For example, bounds on L^{dTi), 
1 < p < 00 , follow from classical Calderon-Zygmund theory, and bounds on layer 
potentials in the domain 12 may be established using bounds on 912; see, for example, 
|DV90| . 

In the case of second-order equations with variable 2-independent coefficients, 
a number of boundedness results have been established. In |KR,09j . Kenig and 
Rule established that in dimension n 1 = 2, layer potentials for operators with 
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real-valued coefficients are bounded on LP{dil) for the domain above a Lipschitz 
graph, and in [R,ul07] this result was extended to boun ded Lipsch itz domains and 
Lipschitz graph domains with arbitrary orientation. In (AAA~*~ iT] , boundedness of 
layer potentials on L^{dfl) was established in arbitrary dimensions, in the domain 
above a Lipschitz graph, for coefficients that are real-valued and symmetric. A 
stability result was also established; that is, if solutions and layer potentials for some 
operator Lq have certain boundedness and invertibility properties on L^(9r2), then 
so do layer potentials for any operator Li whose coefficients Ai are ^-independent 
and near (in L°“) to those of Lq. (This result required a local Holder continuity 
estimate for solutions to Lqu = 0; this estimate is always valid if Aq is real-valued 
but may not be valid for complex Aq.) 

More generally, in [Rosl3] Rosen showed that layer potentials are always bounded 
on L^(9H), for H the domain above a Lipschitz graph, provided that the coefficients 
of the associated operator are t-independent, and also that solutions to Lu = 0 are 
continuous and satisfy the local bound 



whenever Lu = 0 in B{X,r). (The local Holder continuity requirement, used in 
[AAA + n] and in many other papers, is a stronger requirement than this local 
bound. The local boundedness estimate is necessary for Rosen’s construction of 
the fundamental solution E^{X,Y), and thus for the formulas (11.31) and (11.41) to 
be meaningful; he also showed that, even without local boundedness, the double 
and single layer potentials may be continued analytically to bounded operators for 
t-independent coefficients A.) Rosen’s results built on an alternative approach to 
boundary-value problems involving semigroups [AAH081 I A AM 10) : essentially he 
established that layer potentials are equal to certain operators studied in [A AM 10) , 
and thus the boundedness results therein apply. The results of [AAMlOl IRosl3] 
extend to the case of elliptic systems. 

In the case of two dimensions, or of smooth coefficients, standard Calderon- 
Zygmund theory allows for straightforward generalization of bounds to 
bounds, 1 < p < oo. In the case of rough coefficients in higher dimensions, 
new arguments are necessary to bound the layer potentials 0131 and (El on 
L^(9f2). Some such arguments are presented in various papers, in particular in 
[HKMPI 5bl IHMM1 31 IHM!^ . 

In the case of scalar equations, Rosen’s boundedness result was later es¬ 
tablished another way, without semigroups, by Grau de la Herran and Hofmann 
in )GdlHH| . As in I A A A~*~ 1 T] . they required that solutions to Lu = 0 be locally 
Holder continuous, and in particular locally bounded. In this paper we will closely 
follow their approach. We will need to confront a number of additional difficulties 
that arise in the case of higher-order equations. However, one significant advantage 
of the higher-order setting is that local Holder continuity is automatic in the case 
of operators of very high order, and there are established techniques to generalize 
to operators of low or moderate order (see |AHMT0T1 IBarI4) or Section [TT]) ; thus, 
our Theorem [L6l is valid without any assumptions on local boundedness or Holder 
continuity of solutions. 

1.4. Layer potentials for higher order operators: known approaches and 
new ideas. Turning to the history of higher order problems, we recall that an 
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interesting first step lies in even defining layer potentials in the higher order case. In 
particular, the prototypical higher order operator, the bilaplacian A^, can be viewed 
in two ways: either as an operator in the divergence form (EH), A2 = J2j,k9jkidjk), 
or as a composition of two second order operators (two copies of —A). Many 
papers have used potentials based on a formulation of higher order operators as 
compositions; see [DKV 86 i IVerQOi [PV92[ IBM13a) . A somewhat different approach 
is necessary if we view A^ as a divergence form operator, or if we seek to generalize 
to other such operators. 

We begin with Neumann boundary values. Notice that the Neumann boundary 
values V • A*VE^ of the fundamental solution appear in the definition of the 
second order double layer potential. In fact, layer potentials are deeply connected 
to Dirichlet and Neumann boundary values of solutions in other ways; for example, 
if Lu = 0 in n for some second-order operator L, then u satisfies the Green’s 
formula 


(1.9) 


In w = -V^{u\g^) + S^{v ■ AVu). 


That is, such a function u is equal to the double layer potential of its Dirichlet 
boundary values added to the single layer potential of its Neumann boundary values. 

The formulation of Neumann boundary data for higher-order equations is an 
interesting and very difficult question in its own right. It has often been based 
on an integration by parts: for sufficiently nice domains D, operators L, and test 
functions w and v, 


( 1 . 10 ) 


w Lv = 


|a| = |/3|=m 


/ 

Jn 


d°‘w d^v + 


m — 1 

E 

3=0 


/. 


diw Byv da 


for some functions B^v, where dl is the jth derivative in the normal direction. (If 
desired, exact formulas for the functions B-^v in terms of the higher derivatives of v 
may be computed; if L = A^ then formulas for B'^v may be found in |CG85[ fVerOS] 
or in Section [2.2.II below, and in the case of general constant-coefficient operators 
an explicit formula may be found in |MM13bl Proposition 4.3].) 

It is very natural to regard the array as the Dirichlet boundary values 

of w. Then the array of functions {B^v}'JIEq^ may be regarded as the Neumann 
boundary values of v. The Neumann problem for the biharmonic function A^, with 
this formulation of boundary data, was studied in [GG851 IVerOSI ISheOTbl IMM13a) . 
The Neumann problem for more general constant-coefficient operators was studied 
in [VerlOl IMM13b) , and for some classes of variable coefficient operators in |Agr07| . 
We remark that a given operator L may be associated to more than one coefficient 
matrix A, and that each choice of A will give rise to different boundary operators 
B^. We will provide more details and a specific example of these different boundary 
operators (for L = A^) in Section [2.2.11 several of them are physically relevant (in 
different contexts!) and some even lead to ill-posed problems. 
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Going further, if Lu = 0 in SI and {X, Y) is the fundamental solution to L* 
(so that LyE^ then 


In 


(X)u(X)= [ uL*E^^{-,X) 
Jn 


m 

m — 1 


= ^ f diuBfE^^{yX)da-( B 
j=o Jan ' 


^udlE^* {■ ,X) da. 


This naturally suggests the two multiple layer potentials 

m — 1 

IdQ 


m — 1 p 

25 n/W=E/ BfE^‘{-,X)f,da, 

JdQ. 


3=0 
m — 1 


S^giX)=J 2 diE^\-,X)g,da. 
Ian 


3=0 


Notice that in the higher-order case, layer potentials take as input an array of 
several functions. Also, this formulation of layer potentials generalizes the Green’s 
formula (HH). Layer potentials constructed in this way, from an integration by parts 
against the fundamental solution, have been used in [GG831IGG851 IVerOSl IMM13a] 
to study the biharmonic operator (and in particular the associated Neumann 
problem), and in |Agm57[ IMM13bj to study more general constant-coefficient op¬ 
erators; therein boundedness results for layer potentials have been established. 
(Some boundedness results for inputs in fractional smoothness spaces were estab¬ 
lished in |MM13b) .l 

Our formulation of the Neumann boundary values of a solution, and thus layer 
potentials, is different. Specifically, observe that the different terms dlu exhibit 
different degrees of smoothness; if G for example, and dfl is suf¬ 

ficiently smooth, then we expect dlu to lie in the Sobolev space lT^_i_j(9fl) of 
functions with gradients of order m — 1 — j. Furthermore, if G L^(9f2), then 
we generally expect the Neumann boundary terms B'^u to lie in negative smooth¬ 
ness spaces (specifically, we expect B-^u G fFj+i_^(9n), and so only B:^_^u lies 
in L^(9n)). See Section [2 .2.1 1 for an example. 

This is somewhat problematic in the case of Lipschitz and other non-smooth 
domains, as higher smoothness spaces and negative smoothness spaces are difficult 
to formulate. Furthermore, dealing with mixed orders of smoothness is difficult 
even in smooth domains. To avoid these difficulties, we will prefer to regard 
rather than as the Dirichlet boundary values of u; this will allow us to 

formulate a similarly homogeneous notion of Neumann boundary data. The latter 
has the advantage of working with elements of the same degree of smoothness and 
being naturally adaptable to our general context. However, explicit formulas for 
Neumann boundary data can only rarely be obtained; we treat the entire package 
of Neumann data as a linear functional on a suitable Sobolev space. See a detailed 
discussion and an example in Section [2. 2. II We have also formulated layer potentials 
based on this notion of boundary data; see Section 12.41 Our potentials thus take 
as input arrays of functions in homogeneous spaces; notice the norms on the 
right-hand sides of the bounds (11.711 and p.8l) . 
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A Green’s formula involving homogeneous boundary data has been used in 
[PV95bl IVer96] . However, this Green’s formula was formulated in terms of deriva¬ 
tives of order 2m — 1, and as such does not lend itself to formulation of Neumann 
boundary data or the natural division into double and single layer potentials. Fur¬ 
thermore, their construction used some delicate integrations by parts not available 
in the variable coefficient case, and so our formulation of layer potentials is of 
necessity somewhat different and more abstract. 

1.5. Our method and outline of the paper. The remainder of this paper will 
be devoted to a proof of Theorem ll.61 Specifically, we will define our terminology in 
Section[2] We will provide a few preliminary arguments, mainly involving the theory 
of solutions to higher-order equations, in Section [3l We will show that the bounds 
uni) and (11.81) follow from more convenient bounds (specifically, bounds involving 
derivatives in the t-direction only) in Section |4l we will also define new operators 
of and of that are somewhat easier to work with. The proof of Theorem 11.61 
will make extensive use of T1 and Tb theorems; we will state the theorems we 
will need (taken from [CJ87] and |GdlHH] ) in Section [S] The remaining sections 
of the paper will be devoted to showing that Of and Of satisfy the conditions 
of Theorems 15.21 and 15.111 and thus satisfy appropriate estimates; a more detailed 
outline of Sections IMTTl is provided in Section [5j 

Acknowledgements. We would like to thank the American Institute of Math¬ 
ematics for hosting the SQuaRE workshop on “Singular integral operators and 
solvability of boundary problems for elliptic equations with rough coefficients,” at 
which many of the results and techniques of this paper were discussed. 

2. Definitions 

Throughout we work with a divergence-form elliptic operator L of order 2m 
acting on functions defined in R"+^. 

We will reserve the letters a, /3, 7 , ( and ^ to denote multiindices in If 

C = (Cii C 2 , • ■ •, Cn-i-i) is a multiindex, then we define |(C|, 9^ and C! in the usual ways, 

as Id = Cl + C 2 H-f Cn-Li, = di\dil ■ ■ ■ and C! = Cd C 2 ! • • • Cn-hd- K C 

and C are two multiindices, then we say that C ^ C if d ^ Ci for all 1 < i < n -I- 1 , 
and we say that C < C if in addition the strict inequality d < Ci holds for at least 
one such i. 

We will routinely deal with arrays F = (fy;) of numbers or functions indexed by 
multiindices C with |C| = k for some k. In particular, if is a function with weak 
derivatives of order up to k, then we view as such an array. 

The inner product of two such arrays of numbers F and G is given by 

{F,G)= ^ nCc. 

ICI=fc 

If F and G are two arrays of functions defined in an open set H or on its 
boundary, then the inner product of F and G is given by 

= E / or (F,G)g^= E / 

|C|=fe“'^ IC|=fc“'^o 

where cr denotes surface measure. 
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If G is an array of functions defined in n and indexed by multiindices a with 
I a I = k, then div^ G is the distribution given by 

(2.1) (v.,div,G)^ = (-l)™(vV,G)n 

for all smooth test functions (p supported in fl. In particular, if the right-hand side 
is zero for all such (p then we say that divfc G = 0. 

We let Cfc be the unit vector in in the /cth direction; notice that Ck is 

a multiindex with \ek\ = 1- We let be the “unit array” corresponding to the 
multiindex C; thus, (e^, j') = Fc^. We will often distinguish the n -|- 1th direction; 
we let 7 _l = (m — l)e„+i = (0 ,..., 0, m — 1) and a_L = men+i, and let the array 
ej_ denote either or Which is meant should be clear from context. 

We let LP{U) and L°°{U) denote the standard Lebesgue spaces with respect to 
either Lebesgue measure (if t/ is a domain) or surface measure (if C/ is a subset of 
the boundary of a domain). We denote the homogeneous Sobolev space W^{U) by 

Wl{U) = {u : G LP{U)} 

with the norm = ||V^u||iP([/). (Elements of Wj^{U) are then defined only 

up to adding polynomials of degree fc—1.) We say that u G or u G 

if u G L^(E) or u G Wj!{V) for every bounded set V with E C ?7. In particular, if 
[/ is a set and U is its closure, then functions in L^^JU) are required to be locally 
integrable even near the boundary dU] if U is open this is not true of LP^^{U). 

If /r is a measure and E is a /i-measurable set, with p{E) < oo, we let f dp = 
Jl(e) Ie f E C is a set, we let 1e denote the characteristic function 

of E] in particular, we will let 1± denote the characteristic function of the half-space 
If / is a function defined on E, we will often let l^/ denote the extension 
of / to R"+^ by zero. 

If Q C K" is a cube, we let i{Q) be its side-length. We let rQ be the concentric 
cube of side-length r£{Q). We will make frequent use of “dyadic annuli” defined as 
follows. We let 

(2.2) AoiQ) = 2Q, AjiQ) = \ VQ for all j > I. 

If f > 0, let 

i+i 

(2.3) A,^,{Q)= U A,{Q) 

l=i-i 

where A(,{Q) = 0 whenever £ < 0. 

Throughout the paper we will work mainly in the domain = {(a:,t) : a: G 

R",t > 0}. We will also need to consider R"’*'^ = {(a;,t) : a: G R",t < 0}. We will 
often identify R’" with aR±+\ 

If (/? is a function defined on an open subset of R"+^, we will let V||(/? denote the 
gradient only in the first n variables; we will also use V||/ to denote the gradient 
of a function / defined on R" = We will view as an array of functions 

indexed by multiindices C G with |C| = k and Cn+i = 0; equivalently we may 
view V|(/3 as an array of functions indexed by multiindices ^ G N". 

Similarly, if F is an array of functions indexed by multiindices C G N" with 
|(C| = k, we will define div^ y F formally as div^ || F = X)|q=/c ceN^ weak 

definition is precisely analogous to the definition of divfe F. 
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2.1. Elliptic operators. Let A = (Aa^j be measurable coefficients defined on 

indexed by multtiindices a, /3 with |a| = |/3| = m. If E is an array, then AF 
is the array given by 

(AF)„= ^ 

|/3|=m 

Throughout we consider coefficients that satisfy the Carding inequality 

(2.4) Re > A|| for all <p G 

and the bound 

(2.5) II A||j^oo(][{n + l) < A 
for some A > A > 0. 

We let L be the 2TOth-order divergence-form operator associated with A. That 
is, we say that Lu = div^ E in fl in the weak sense if, for every (p smooth and 
compactly supported in fl, we have that 

(2.6) (V-^,AV-ir)^ = (V-^,E)^, 

that is, we have that 

(2.7) ^ f d^^A^pd^u= Y, 

In particular, if the left-hand side is zero for all such ip then we say that Lu = 0. 

We let A* be the adjoint matrix, that is, We let L* be the associated 

elliptic operator. 

In this paper we will focus exclusively on operators L that are t-independent, 
that is, whose coefficients satisfy formula dm. 

Throughout the paper we will let C denote a constant whose value may change 
from line to line, but which depends only on the dimension n -I- I, the ellipticity 
constants A and A in the bounds (j2.4p and (12.5p . and the order 2m of our elliptic 
operators. Any other dependencies will be indicated explicitly. We say that A k, B 
if, for some such constant C, A < CB and B < CA. 

2.2. Dirichlet and Neumann boundary data. Our goal in the present paper is 
to bound the double and single layer potentials; in future work we hope to use the 
results of this paper to solve the Dirichlet and Neumann boundary value problems. 
Thus, in this section, we will define higher-order Dirichlet and Neumann boundary 
data. 

We define higher-order Dirichlet boundary data as follows. Suppose that V™E G 
or V'^E G Then d'^u G fEi^^(R^) for any 7 with 

I 7 I = TO — I. We define Tcm-i u as the array given by 

( 2 . 8 ) (Tr^-i u)^ = Trd'^u for all I 7 I = to — I 

where Tr denotes the standard trace operator on Sobolev spaces. In some cases we 
will use Tr^_i u and Tr~_]^ u to specify that the trace is taken in or 

Notice that if 'V'^u is locally integrable in all of R"+^, then u = u. 

With some care we may define boundary values of certain higher-order deriva¬ 
tives. If M G RA that Tcm-i u G kki/oc(®")> then for each /3 with 

Pn+i < |/3| = m, we define 

(2.9) (T^m.l u)f 3 = dxj Tr(9^“®^M for all j with 1 < j < n and Pj > 0. 
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Note the requirement that j ^ n + 1. This is well-defined; that is, if /3j > 0 and 
/3fc > 0 for some j ^ n + 1 and k ^ n -I- 1, then it does not matter whether we 
choose Xj or X}^ as our distinguished representative. 

We define higher-order Neumann boundary data as follows. Let IT,^ q(IR 5;^^) 
be the closure in of the set of all smooth functions supported in 

Then 

^™.o(R±+') = {^ e : Tr±_i v = 0}. 

Observe that if u G IT^ u G and Lrt = 0 in then 

(V’"i;,AV™w)k.+i =0 . 

Thus, if (/? G then the inner product ylV'"'u)gn+i depends only 

on Tri^_i if Tr^-i ^ — "^m-i ^ —V™??, plV’”M)gn+i = 0. We define 

the Neumann boundary values u by 

(2.10) (Tr±_i^,M± for all G lT,?,(Rr')- 

u is then a linear operator on the space of traces of W^(R£^^)-functions. 

2.2.1. Historical remarks and context. We now provide some further discussion and 
history of Dirichlet and Neumann boundary data. 

We remark that we have three ways to refer to derivatives of order m at R" = 
9R![.^^, namely V™m(x, 0), u{x), and V™u(x, 0). All three are arrays indexed 

by multiindices /3 with |/3| = m. To give the reader some intuition, let us discuss the 
simplest case, that of smooth functions in R^ = {{x,t) : x G R, t > 0} with m = 2. 
In this case, formally, Tri u is the array containing dxu{x, 0) and dtu{x, 0). As for 
the higher order traces, V'^u{x, 0) = dl^u{x, 0), Tr„ | u{x) is the array d^^u) 

containing = dx TidxU and dxtU = dx Trdtu on 9R^, while V'^u{x, 0) is the 
array {d^xU, d^^u, dffu) of all second derivatives. 

The reader should compare our choice of representation of the Dirichlet data 
Tri u = (dxu(x, 0), dtu(x, 0)) to the traditional choice (u(x,0),dtu(x,0)). This is, 
of course, a question of representation and what matters is the function spaces 
for the data. Working with Tri in place of (u(x,0),dtu(x,0)) brings considerable 
advantage and clarity mainly because both parts of the trace array belong to a 
function space of the same level of smoothness. For example, Tri is a vector with 
both components in L^(R) when (u(x,0),dtu(x,0)) lies in IF^’^(R) x L^(R). This 
makes things much clearer when dealing with divergence form operators of arbitrary 
higher order. This also allows us to properly define Neumann data. 

Recall that the latter is quite tricky even for the simple case of the bilaplacian 
and that some choices can make the Neumann problem ill-posed. Let us discuss 
this in some detail, starting with the bilaplacian on a Lipschitz domain D C R". 
We will translate to the half-space below. 

The Neumann boundary values of a solution are traditionally given by an integra¬ 
tion by parts (formula (ll.lOp l or less explicitly as an inner product (formula (12.101) 1. 
In the case of the biharmonic equation, Neumann boundary values also have ap¬ 
plications in the theory of elasticity. Recall that the principal physical motivation 
for the inhomogeneous biharmonic equation A^it = h is that it describes the equi¬ 
librium position of a thin elastic plate subject to a vertical force h. The Dirichlet 
problem = /, = g describes an elastic plate whose edges are clamped, 
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that is, held at a fixed position in a fixed orientation. The Neumann problem, 
on the other hand, corresponds to the case of a free boundary. Guido Sweers has 
written an excellent short paper [Swe09] discussing the boundary conditions that 
correspond to these and other physical situations. 

More precisely, if a thin two-dimensional plate is subject to a force h and the 
edges are free to move, then its displacement u satisfies the boundary value problem 

A^u = h in n, 
pAu -b (1 — p)d1u = 0 on dVl, 
d^Au -b (1 — p)drTuU = 0 on 90. 

Here p is a physical constant, called the Poisson ratio, and v and r are the unit 
outward normal and unit tangent vectors to the boundary. This formulation goes 
back to Kirchoff and is well known in the theory of elasticity; see, for example. 
Section 3.1 and Chapter 8 of the classic engineering text [Nad63) . We remark that 
by |Nad63l formula (8-10)], 

dvAu -b (1 - p)drTv'U = d^Au -b (1 - p)dr {dvru) . 

This suggests the following homogeneous boundary value problem in a Lipschitz 
domain fl C of arbitrary dimension. We say that the L^-Neumann problem 

is well-posed if there exists a constant C > 0 such that, for every /o S U’{d^) and 
Aq e lT^j(9fl), there exists a function u such that 


( 2 . 11 ) 


A^u = 0 

MpU := pAu -b (1 - p)dlu = fo 
Kpu := dyAu -b (1 - (9 


^ ||fV(V^M)||iP(ao) < (911/011^^(30) + (9||Ao||wfj(3n) 


in f2, 
on dQ, 

) = Aq on 90, 
on 90. 


Here = ViCj — Vjti is a vector orthogonal to the outward normal v and lying in 
the a;iaij-plane, and denotes the nontangential maximal function common 

in the theory of elliptic boundary value problems. We apply the convention that 
the repeated indices j and k are summed from 1 to n -b 1. 

The boundary operators Mp and Kp, derived from the theory of elasticity, are 
the same as the Neumann boundary operators discussed in Section ri.4l Specifically, 
for any p £ R, the equation 


( 2 . 12 ) 


’ A'^v = / {p Aw Av + {1 — p) djkw djkv) + 


I ^ 


w KpV — dyW MpV da 


is valid for arbitrary smooth functions. Comparing to formula (II.IOL we see that 
= Kp and = —Mp, where A = Ap is an appropriate choice of coefficients 
associated with L = A^. 

Observe that, contrary to the Laplacian or more general second order operators, 
there is a family of relevant Neumann data for the biharmonic equation. Moreover, 
different values (or, rather, ranges) of p correspond to different natural physical 
situations. We refer the reader to [Ver05| for a detailed discussion. 

Recall that our formulation of Neumann data is somewhat different; we use the 
array MJ u of formula (12.101) rather than the functions B-^u of formula (11.101) . As 
an example, in the case L = and H = we will provide an explicit formula 
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for one representative of M^u. In the half-space, 

Kpu = -dn+iAu - (1 - 

We are now summing from j = 1 to n. If Au is harmonic in then dn+iAu = 

djRj{Au), where Rj denotes the jth Riesz transform. Thus, if (p : R" i—>■ C" is any 
divergence-free vector field, then 

(2.13) 

/ wA^u= / {p Aw Au + {1 — p) djkW djku) 

+ / dn+iwMpU + djw{Rj{Au) + {l-p)d'^,^_^_^)U + ipj)da 

from which we may recover an explicit formula for u. 

We comment on several aspects of this formula. First, observe that we still have 
a family of Neumann boundary data indexed by the parameter p. Next, observe 
that we did use the fact that A^u = 0; this formula, unlike formula (12.121) . is not 
valid for arbitrary smooth functions. Furthermore, observe the presence of the 
vector field 0 in m; our explicit representation gives a natural normalization 
(p = t), but for more general operators the divergence-free vector field cannot be 
neglected. Finally, observe that our formula for u is not a local one: it involves 
the Riesz transforms of derivatives of u rather than simply linear combinations. 

However, notice one significant advantage of our formulation (12.131) over the 
operators Kp and Mp. The term MpU involves second derivatives of u, while the 
term KpU involves third derivatives; we have expressed all components of u 
using the second derivatives of u. As discussed in Section [T4l this means that we 
expect the different components of the Neumann boundary data to lie in a single 
smoothness space; furthermore, using boundedness of the Riesz transform, we may 
control ||M(^M||ip(Rn) by ||V2u||ip(Rn), for 1 < p < oo. 

Our formulation of Neumann boundary data for general operators will display 
most of these issues. The existence of a family of Neumann data may be eliminated 
by specifying the matrix of coefficients A in formula (12.101) , but our formulation of 
M+ u does require that m be a solution, is well-defined only up to adding divergence- 
free terms, and need not have a local representation. Indeed, at this point, the 
estimate ||M)4 < CH V™u||iP(Rn), while apparently plausible, is still only 

a conjecture. 

We now discuss the history of the L^-Neumann problem (12.111) . In [CG85] . Cohen 
and Gosselin showed that this problem was well-posed in domains contained in 
for for 1 < p < oo, provided in addition that p = —1. In [VerOSj . Verchota 
investigated the Neumann problem (12.111) in full generality. He considered Lipschitz 
domains with compact, connected boundary contained in n -|- 1 > 2. He 

showed that if — 1/n < p < 1, then the Neumann problem is well-posed provided 
2 — e < p < 2 + e. That is, the solutions exist, satisfy the desired estimates, 
and are unique either modulo functions of an appropriate class, or (in the case 
where H is unbounded) when subject to an appropriate growth condition. See 
[VerOSl Theorems 13.2 and 15.4]. The Neumann problem is ill-posed for p > 1 
and p < 1/n; see [Ver051 Section 21]. More recently, in [She07a) . Shen improved 
upon Verchota’s results by extending the range on p (in bounded simply connected 
Lipschitz domains) to 2n/(n -f 2) — e < p < 2 -f e if n -I- 1 > 4, and 1 < p < 2 -|- e 
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ifn + l = 2orn + l = 3. All of the aforementioned results rely on the method of 
layer potentials. Finally, in [MM13bl Section 6.5], I. Mitrea and M. Mitrea showed 
that if n C is a simply connected domain whose unit outward normal v lies 

in VMO{dVl) (for example, if 11 is a domain), then the acceptable range of p is 
1 < p < oo; this may be seen as a generalization of the result of Cohen and Gosselin 
to higher dimensions, to other values of p, and to slightly rougher domains. The 
question of the sharp range of p for which the L^-Neumann problem is well-posed 
in a Lipschitz domain is still open. 

It turns out that extending well-posedness results for the Neumann problem 
beyond the case of the bilaplacian is an excruciatingly difficult problem. Even 
defining Neumann boundary values for more general operators is a difficult problem; 
compare the traditional definition of Neumann boundary data of (|1.10p and our 
formulation (I2.10p . While some progress has been made (see |Agr07|lMM13bl[Baj] i. 
at present there are no well-posedness results for the Neumann problem with 
boundary data. 

Recall that |Agr07[ IMM13b] have investigated the Neumann problem for bound¬ 
ary data formulated as in formula (11.1011 . Specifically, Agranovich has established 
some well-posedness results for the inhomogeneous problem Lu = h with homo¬ 
geneous Neumann boundary data, and has provided some brief discussion of the 
conditions needed to resolve the Neumann problem with inhomogeneous boundary 
data; see |Agr07[ Section 5.2]. The book |MM13b) considers the case of constant- 
coefficient operators at length; therein they establish well-posedness results for the 
Neumann problem, with boundary data in certain fractional smoothness spaces, for 
elliptic constant-coefficient operators. 

2.3. The Newton potential and the fundamental solution. The main pur¬ 
pose of the present paper is to define and bound the double and single layer poten¬ 
tials for higher-order elliptic operators of the form specified in Section \TJ\ Recall 
from formulas (O and dn that the second-order layer potentials are built from 
the second-order fundamental solution. 

The main result of the paper [Bar 14] was a construction of the fundamental 
solution in the case of higher-order operators. was constructed as (an 
order-TO antiderivative of) the kernel to the operator 11^, the Newton potential 
for L, defined as follows. For any H G L^(M"“''^), by the Lax-Milgram lemma there 
is a unique function Il^H in (]&"■+^) that satisfies 

(2.14) AV™n^iT)K„+i = (V™(^, 

for all (fi € kF^(R"+^). The Newton potential is a bounded operator on T^(IR"+^) 
and satisfies the bound 

(2.15) < C'||JT|U2(r„+i). 

We will need two additional properties of the Newton potential. First, we will need 
the symmetry relation 

(2.16) (G, V’-n^iT)Rn+i = (V™n^*G, iT)r+i 

for all H G and all G G L^(K"). Second, we will need the identity 

(2.17) V"^n^(AV"^E) = 

for all E G 1F^(M"“''^); this identity follows by uniqueness of the Newton potential 
as the solution operator in formula (I2.14|) . 
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We remark that this Newton potential 11^ is somewhat different in smoothness 
from the traditional Newton potential. This potential, which we will denote N^, is 
the unique solution of the equation Lu = f, u G or, more precisely, of 

(2.18) = (tp, /)«„+! 

for all (f G The input / should thus be taken in the dual 

space to One can write / = div^ for some jj G The 

above formulas then become Lu = divm H, u G or, more precisely, 

(2.19) AV™iV^ div^ H)un+i = iT)R„+i. 

In other words, = N^^divm- We shall be working exclusively with n-^, but 
perhaps this analogy is useful to keep in mind. 

The main result of [Barl4| may be stated as follows. 


Theorem 2.20 ( |Barl4l Theorem 62 and Lemma 69]). Let L he an operator of 
order 2m that satisfies the bounds (EH) and ( 1 ^ . Then there exists a function 
E^{X,Y) with the following properties. 

Let q and s be two integers that satisfy q + s < n + 1 and the bounds 0 < q < 
min(TO, {n + l)/ 2 ), 0 < s < min( 77 i, (n + l)/ 2 ). 

Then we have the symmetry property 


(2.21) dj^dlE^{X,Y) = dj^dlE^" {Y,X) 

as locally functions, for all multiindices f with \(^\ = m — q and j^j = m — s. 

There is some e > 0 such that if Xq, Yq G if 0 < dr < R < \Xo — To|/3, 

and if q < (n + l)/2 then 

(2.22) [ [ |V'^"®V™“'^L;^(X,y)|2dXdy < 

JB{YQ,r) JB(Xq,R) 

If q = (n + l)/2 then we instead have the bound 

(2.23) [ [ \Xx-’'XY~'^E^{X,Y)\^dXdY <CiS)r^'^R^’’ 

JB{Yo,r) Jb(Xo,R} 

for all 6 > 0 and some constant C{d) depending on <5. 

Furthermore, if \a\ = m then 




(2.24) 



d^d^E^{X, Y)Fp{Y)dY 


for almost every X ^ suppiT, and for all H G L^(M"“'"^) whose support is not all 
of R"+i. 

Finally, if E^ is any other function that satisfies the hounds (12.2211 . (12.231) and 
formula (I2.24p . then 

(2.25) E^{X,Y) = E^{X,Y) + E fdX)Y<+ E 

\C\<m—{n+l)/2 |{|<m —(n+l)/2 

+ E 

ICI = l?l=”i-("+l)/2 


for some functions and gc^ and some constants C(^^. Thus, Y) 

is a well-defined, locally function provided q and s satisfy the conditions specified 
above. 
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Note that formula p.24|l and the definition of assures that E is indeed 
analogous to the traditional fundamental solution, which, roughly speaking, solves 
LE = 5. That is, E is formally the kernel of the potential N defined above. 

We record one further property of the fundamental solution for t-independent 
operators. By the uniqueness property for the fundamental solution, if A is t- 
independent, then we have that 

di,t9isE^ix,t,y,s) = di^td^^^E^{x,t + r,y,s + r) 

for almost every x € R", y G R" and almost every t, s, r G R, and all multiindices 
^ as in formula (j2.21ll . In particular, for such C and ^ we have that 

(2-26) di^tdl^,dtE^{xA,y,s) = -di^tdl^,d,E^{x,t,y,s). 

Remark 2.27. We comment on the additional terms in formula (I2.25|) . Notice 
that E^ is defined essentially by the relation (12.241) . But this relation involves only 
derivatives of order 2m; in other words, it is only V^V^E^{X,Y) that is well- 
defined. The lower-order derivatives are defined only up to adding polynomials. 
(The dx derivative is included in formula (|2.24l) because U^H G IT^(R"+^), and 
so n^H is also defined only up to adding polynomials.) 

If q and s are small enough, then there is a unique normalization of the deriva¬ 
tives y) that satisfies the bound (12.2211 or (12.2,41) : in [Bar 14] 

this normalization was found using the Gagliardo-Nirenberg-Sobolev inequality. 
However, if 2m > n -|- 1 then E^ itself (and possibly some of its derivatives) are 
still not well-defined. The extra terms on the right-hand side of formula (12.251) 
are precisely the terms compatible with the requirement T) = 

y) for q, s small enough. 

Consequently, throughout this paper we will be careful to use only derivatives 
of E^ of sufficiently high order; in fact, we will use only derivatives of the form 
dxdyE^iX, Y) for |^| > m — 1, |C| > m — 1 and |5| -I- |C| > 2m — 1. 

In some very special cases, there are natural normalization conditions for the fun¬ 
damental solution even if 2m > n -|- 1; for example, if L = (—A)™ and n -|- 1 < 2m 
is even, then we may take E^{X,Y) = Cm,n\X — y|2m-(™+i) log|A — y|. Notice 
the presence of logarithmic growth in the fundamental solution. However, if we 
take 2m — (n -I- 1) -I- 1 derivatives (in either X or y), then the logarithm vanishes; 
this is the lowest order of derivative that Theorem 12.201 guarantees is well-defined. 

2.4. The double and single layer potentials. In this paper we seek to formulate 
a notion of layer potentials for higher-order elliptic operators of the form specified 
in Section o The goal of this paper is to produce bounds on layer potentials in 
the domain H = R!(.^^; thus, we will define boundary values and layer potentials 
only for the half-spaces. 

We begin by recalling the second-order Green’s formula (HH). To generalize this 
formula to higher order, notice that, for any function u G iy^(R"^^), 

1+u = (1+u - n^(l+AV’"w)) -b n^(l+AV’"u) 

as iy^(R!(.^^)-functions. We claim that the quantity 

(2.28) =-l+F-l-n^(l+AV™y) if/ = Th+_iF 

is well-defined; that is, the right-hand side depends only on E. We will define 

S^g in such a way that 5^(M)^u) = n^(l+AV'"u) as iy^(R’"+i)-functions for 
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any u S with Lm = 0 in this then yields the higher-order Green’s 

formula 

(2.29) Ijj.+i V"*!! = u) + u). 

We will also find formulas (12.311) and (12.351) for the double and single layer potentials 
in terms of the fundamental solution; these formulas will also parallel formulas (131) 
and (jl.4L 

We now establish our claim that, if / = Tfm-i F for some F G then 

V^f is a well-defined element of and It suffices to show 

that, if = 0, then the right-hand side of formula ()2.28l) is zero. 

Suppose that Tr^-i F = Tr = 0. It is well known (see, for example, the 

proof of |Eva981 Theorem 5.5.2]) that in this case, lies in the completion in 

Wi^(]R"~'"^) of the set of smooth functions compactly supported in Now, if ip 

is compactly supported in then we may extend p to a, function in all of 

by letting tp = 0 in By density, we have that also extends by zero to 

a function in Wf (IR"+^); thus, we may extend F to a polynomial of degree m — I 
in Without loss of generality we may take this to be the zero polynomial. 

Thus, if F = 0, then l+F G W^(]R"+^). We may apply the identity (12.171) 

to 1+F, and so the right-hand side of formula (12.281) is zero in that is, 

up to adding polynomials of degree to — I. 

We will need two alternative formulations of Notice that we may extend 

F to a function even if Tr^_i F = / ^ 0; then Tr~_i F = / as well. 

Then by formula (I2.17p . 

(2.30) = l_F-n^(l_AV'"F) if / = F. 

By formula (12.241) . if jaj = to, then for almost every x G M" and t > 0, we have 
that 

(2.31) d°^V^f{x,t) = - ^ ^^^^l^E^{x,t,y,s)Ap^{y,s)^^F{y,s)dsdy. 

|/3|=m ®- 
|5|=m 

A corresponding formula, involving an integral over is valid if f < 0. 

We will establish a bound on in terms of the norm of the tangential 

derivative Vy/ of /. In order to use existing theorems concerning boundedness, 
we will want to slightly modify the definition of the double layer potential, by 
defining 

(2.32) F^(Tr™.| F)(x, t) = F)(a;, t) 

for all sufficiently well-behaved functions F. 

We now must define the single layer potential. Let g be a bounded linear operator 
on the space 

W4 ^_i/2 (R”) = {fr+_, F : F G IF^(R(^+')}. 

The operator TgF = {g, Tr^_-^ F)^n is a bounded linear operator on 
We may identify IT^(R"'''^) with a subspace of (F^(R"^^))'J, where q is the number 
of multiindices a of length to, via the map F i—>■ V™F. We may then extend Tg 
to an operator on (F^(R"“''^))‘^. Let G G (F^(R"^^))'^ be the kernel of Tg, so 
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Tg{H) = (G,A)Kn+i for all H € (L2(IR^+1))9, in particular, (G, = 

(g,Tr+_i V5)aR:;;+i for all ip S W^(R++^). Let 

(2.33) S^g = n^(l+G) if (G, = {g, A+_i for all p e W^. 

As in the case of the double layer potential, it is straightforward to establish that 
S^g S and that the value of S^g is independent of the choice of G, 

that is, of the choice of extension of Tg from to The for¬ 
mula = n^(l+AV"*M) follows immediately from the definitions (12.101) 

and (12.331) of Neumann boundary data and the single layer potential. 

Again we will need alternative formulations of layer potentials. First, observe 
that 

(2.34) S^g = n^(l_G) if (G, V'^p)^^_+i = {g, Tr^_i for all p e 

By formula (12.241) . if |a| = m then 

d‘^S^g{x,t) = ^ [ d‘^td^,,E^{x,t,y,s)G0{y,s)dsdy. 


But by the bound (I2.22L p{y,s) = d^fE^{x,t,y,s) is a TF^(R"^^)-function for 
almost every x G M" and t > 0; thus we may write 


(2.35) 


d^S^g{x,t)^ 


E 

l-yl^m —1 


dx,td^y,sE^{x,t,y,Q)g^(y) dy. 


2.5. Function spaces on the boundary. We have now defined and as 
operators on and its dual space, respectively. We wish to extend 

and 5^ to bounded operators on the space L^(]R"). However, notice that acts 
naturally only on traces of gradients; that is, density arguments will only allow us 
to extend to a subspace of L^(]R"). We will define this subspace as follows. 

Definition 2.36. We let WA‘^_i{W^) be the completion of the set 
{Trm_i p : p smooth and compactly supported} 
under the norm. 

We let WA^ | (R") be the completion of the set 

2 ) = {Tr,„j p : p smooth and compactly supported} 
under the norm. 

It is well known that that the space used above is the com¬ 
pletion of 2) under the norm in the Besov space This space is of¬ 

ten called a Whitney-Besov space and has been used in the theory of higher- 
order boundary-value problems; see, for example, [AP981 [Agr07[ IMMS101IMMW 11 1 
IBM13bl lMM13b[ IBar] . The spaces lib4^_;^(R") or WA"^ |(1^") are called Whitney- 
Sobolev spaces; they have also been used extensively in the theory, for example, 
in [Ver901 IPV95al IPV95bl IVer961 LSheOBal ISheOBbl IKS 11) . The goal of this pa¬ 
per is to extend the double and single layer potentials to bounded operators on 
Whitney-Sobolev spaces by establishing boundedness results. 
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Remark 2.37. We remark that is a well-defined operator on the space 
D = {Trm.l <^5 : is smooth and compactly supported}. 

we will extend to WA^ by density. 

We will also extend the single layer potential; in this case we wish to extend 
to all arrays of functions g G L^(K"). It will be convenient to have a dense 
subspace ITl at our disposal on which is known to be well-defined. We claim that 
S^g is well-defined for any g G L^(]R"') that is compactly supported and integrates 
to zero. 

The reasoning is as follows. Recall that is well-defined whenever g is a 

bounded linear operator on 

W4 ^_i/2 (K”) = {Tr„-i $ : G 

Now, suppose that J — 0. Choose some function $ G that is smooth 

up to the boundary. Then 

[ g^i (5^^ - c®) 

dR" Jr" 

for any constant c<i>. Suppose that is supported in R" fl B{{xq,0),R) for some 
xq G M" and some R > 0. Let ft = fl B{{xo,0), R)- It is well known (see, 
for example, |Eva98 ] ) the trace map is bounded from fl Wf{n) to L'^{dfl). 

Thus, 


/ 

JR'- 


9j 


97$ 


< ||57llT^(R")ll^"’’‘b - C$||L2(af2) 

< — C$||L2(n) 




By the Poincare inequality, if we choose c$ correctly then we may control the 
quantity - c<[,||i 2 (o) by ||Va^$|U 2 (o) < ||V™$||^ 2 (Rn+i), and so we see 

that g gives rise to a bounded operator on Thus, for such g, S^g 

is a well-defined element of IT^(K"+^). 


3. Preliminary arguments 

In this section we will establish some basic results that will be useful throughout 
the paper. 

We begin with some bounds on solutions to elliptic equations. Specifically, we 
begin with the following higher-order generalization of the Caccioppoli inequality; 
in its full generality it was proven in [Barl4) . but the j = m case was proven in 
[CamSOj and an intriguing version appears in |AQ00| . 

Lemma 3.1 (The Caccioppoli inequality). Suppose that L is a divergence-form 
elliptic operator associated to coefficients A satisfying the ellipticity conditions dm) 
and (12.51) . Let u G IT^(R(Xo, 2r)) with Lu = 0 m B{Xo,2r). 

Then we have the bound 

-/ \X^u{x,s)\^ dxds < ^-f \V^~^u{x,s)\^ dxds 

J B(X,r) J B(X,2r) 

for any j with 1 < j < m. 
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We may use the following lemma to bound u not in balls of dimension n + 1, but 
on horizontal slices of dimension n; the second-order case of this lemma is known 
and is Proposition 2.1 in |AAA~*~ll] . 


Lemma 3.2. Let t be a constant, and let Q C K" be a cube. 

Suppose that dsu{x, s) satisfies the Caecioppoli-like inequality 

-r \dsu{x,s)\^ dxds < ^ ^ \u{x, s)\^ dx ds 

J B{X,r) J B{X,2r) 

whenever B(X, 2r) C {{x, s) : x € 2Q, t — £{Q) < s < t + £{Q)}. 

Then 

\u{x,t)\‘^ dx </ \u{x,s)\'^ dsdx. 

Jq J2QJt-l{Q) 

In particular, if Lu = 0 m 2Q x {t — £{Q),t + £(Q)), and L is an operator of 
order 2m associated to t-independent coefficients A, then 

rt+t.{Q) 


r c r 

/ dfU{x,t)\'^ dx < / / dgU{x, s)]"^ ds dx 

Jq ~ ^(Q) J2QJt-i{Q) 

for any 0 < j < m and any integer k > 0. 


Proof. Begin by observing that 

[ \Mx,t)\ 

\Jq 


1/2 

^ dx\ < 


Id 


t+e{Q)/2 


u{x, s) ds 


2 s 1/2 

dx I 


t+eiQ)/2 ^ 1/2 

|^t(x, s)p ds dx 


But 


/ u{x,t)-f 
JQ Jt 


t+l{Q)/2 2 

u(x, s) ds 


• r rPQ)/'^ rs 

dx < j d / dru{x,t + r) dr ds 

Jq Jo Jo 

f rPQ)/-^ 

< / / \dru{x,t + r)\^ dr dx. 

Jo Jo 


dx 


Applying the Caccioppoli inequality completes the proof. 


□ 


Throughout this paper we will frequently need to bound the fundamental so¬ 
lution of Theorem 12.201 on horizontal slices. The following estimate follows from 
Lemma [XU the Caccioppoli inequality and the bound (12.221) : we will use it several 
times. Suppose that Q is a cube and that either j > 1 or j > 0 and i{Q) < |s — t\. 
Suppose further that q, s, i and k are nonnegative integers with q < m, s < m and 
q — k < (n -I- l)/2, s — i < (n -|-1)/2. Then for some e > 0, 


(3.3) 


QJA,(Q) 


\VTfi'^V^-^d’^dlE^ix,t,y,.s)\^dydx< 


C 




2-A2{i-s) + l+e) 


where r = k + i — q — s = (m -(/)-!- (to — .s) + k + i — 2m. In applying this formula 
it is always useful to remember formula p.26|l . that is, that we may take vertical 
derivatives in either the s variable or the t variable. 

Now, recall that we seek to bound the double layer potential T>^. Furthermore, 
recall that if / = Trm-iF, then we may write in terms of F; see formu¬ 

las (12.2811 or (|2.31|) . Thus, we will be concerned with extensions of elements of 
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(R"), or of elements of the dense subspace S of Section [231 The following 
lemma provides some basic extensions and some estimates upon such functions. 

Lemma 3.4. Let f = Tr^-i F for some smooth, compactly supported function F. 

Then there is some function FI defined in such that ThCm-i FI — f and 

such that 

(3.5) < C / 

(3.6) sup||V™"iiJ(.,t)||L(Rn) < / \f{x)\^dx, 

t/o Jr" 

(3.7) sup||V-77(.,t)||L(„„)< / |V||/(cr)|2dx, 

Jr" 

(3.8) [ r\V^H{x,t)\ 

Jr^ Jo 

Furthermore, if f = 0 in some cube Q, then = 0 in {{x,t) : dist(a;,K" \ 

Q) > t}, and in particular in {1/2)Q x (0,t'((5)/4). 


'^tdtdx< f \f{x)\'^dx. 

Jr 


Proof. For each 0 < j < m — 1, let fj{x) = di^_^iF{x, 0); observe that up to adding 
polynomials of appropriate degree, fj is determined entirely by / = 0). 

Let ?7 : R" !->• R be smooth, nonnegative, supported in 73(0,1), and satisfy 
/r" V = ^ and x'’ r]{x) dx = 0 for all multiindices C G N" with 1 < |CI < ™ — 1- 
Let ptix) = t~'^r]{x/t). Let 

Hj{x, t) = —F fj * rjtix) = [ fj{x - ty) r]{y) dy, H{x, t) = 

J- J- Jr^ 

By inspection, F[j{x,t) = fj(x), and if0<fc<m — 1 with j k, then 

limt^o dtdlj{x, t) = 0 if j ^ k. Thus Trm_i iJ = /, as desired. 

We may bound FI in terms the functions fj using the Fourier transform in the 
x-variable and Plancherel’s theorem, and we may bound appropriate derivatives 
of fj using the array /. We omit the routine details. 

If / = 0 in Q, then W"^~^~^fj = 0 in Q, and so fj is a polynomial in Q. Thus, 
we may write fj{x — ty) = E| 7 |<m-i-j A(^) some polynomials P-y{x). 

Notice Pq{x) = fj{x). By our moment condition on p, if dist(x,M" \Q) > t, then 

Hj {x,t) = -F [ f (x - ty) T]{y) dy = ^F f Pq (x) r]{y) dy = ^F fj (x). 

J- Jr" J- 7r" J- 

Thus, H(x,t) is equal to a polynomial of degree at most to — 2 in this region, as 
desired. □ 


^ Hj{x,t). 
j=o 


4. Operators to be bounded 

Recall that Theorem 11.61 involves bounding the quantities V'^dn+iS^g and 
V^dn+i'D'^f. In this section we will reduce to the case of the purely vertical 
derivatives; that is, we will show that bounding and for any 

A: > 1, suffices to bound and V'^dn+i'D'^ f. We will also establish 

some notation for these operators. 
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Let A: > 1 be an integer, to be chosen later. Let 

(4.1) Qlg{x)=t^dr+^S^g{x,t). 


Observe that by formula p.35|l . 

(4.2) Qfg{x)= [ d^+^dl,E^{x,t,y,0)g-^{y)dy. 

|7|=m-l 


Notice that by the bound (13.3L if k is large enough and if g G L^(IR."+^), then the 
integral converges for almost every {x,t) € We will elaborate on this point 

in Section [6l 

If f lies in the space ® of Remark 12.371 we let 

(4.3) 


Establishing a bound on 0^ of some sort in terms of the norm of / will allow 
us to extend 0^* to all of WA^ |(R"')- 

Note that 0f, Of implicitly depend on fc > 1. 

We begin by reducing the proof of Theorem 11.61 to establishing bounds on 0f 
and of; the remainder of this paper will be devoted to establishing these bounds. 

Remark 4.4. The conclusion of Theorem II.61 is a bound in the whole space 
for notational convenience, we will establish a bound only in the upper half-space 
R"’*’^ and note that the corresponding bound in follows by careful argument 

involving the change of variables (x,t) i—>■ (x, —t). 


Lemma 4.5. Let f G D and g G fh, where (D and 91 are as in Reniark \2.37\ If 
k > 1, then we have the bounds 



'\/"^dtS^g{x,t)ftdxdt < C 



V^dtV^j{X)ftdxdt<C 



\^t9ix)f\dxdt, 
\Of f{x)f' - dx dt. 


Proof. We follow the proof of the similar formula (5.5) in [AAA~*~lT] . Let u = S'^g 
or It = T>^f, and define 


^jf)— f |V™cl/zt(a;,t)p dx, Vj(t) = f u{x,t)f dx. 

To prove the lemma we need only establish the bound 

pOO poo 

/ tUi{t)dt<C t^^-^Vk{t)dt. 

Jo Jo 


By Lemma 13.21 and the Caccioppoli inequality, if J > — m then 

r2t 


C f 

and thus we may easily show that 


pOO pOO 

/ t“'=-it/^+fc(t)dt < / t^'^-^Vkit) dt. 

Jo Jo 
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Thus, we need only show that 

poo poo 

/ tUi{t)dt<C Um+kit) dt. 

Jo Jo 


Observe that G L 
we have that if j > 0, then 


2(mn+l 


^ ). By Lemma [3.21 and the Caccioppoli inequality, 


C 


uAt) < 

Suppose j > 0. Then ifO<e<S'<oo, we have that 


pS pS pS pS 

J f^-^Uj{t)dt = J t‘^^-^Uj{S)dt- J j u'{s)dsdt 


+ s^^\U'{s)\ds. 


S" 2j 

Observe that |C/'(s)| < 2yJUj{s) Uj+i{s) < ^Uj{s) + sUj+i{s). Thus, 

U,{t) dt < ^11 ^ W ds 

+ Yjj^ C/j+i(s)ds. 

Rearranging terms, we have that if j > 1 then 

nS 


U,{t)dt < u^^^^s)ds. 

Taking the limit as e —>■ O'*' and S' —>■ oo, we have that if j > 0 then 

pOO pOO 

/ t^j-i dt<c Uj+i{s) ds. 

Jo Jo 

Iterating, we see that 

poo poo 

/ tU,{t) dt < C{k) / t2m+2fe-l U^+k{t) , 

Jo Jo 


I dt 


as desired. 


□ 


Remark 4.6. In Sections ISHTTl we will bound the operators and Of for some 
values of A:. In particular, we will make many arguments that are only valid for 
k large enough; we will not make any arguments that are only valid for k small 
enough, and thus there will be some k large enough that all our arguments are 
valid. 


5. A VECTOR-VALUED T{b) THEOREM 

Our goal now is to produce square-function estimates for the operators Of and 
of. In this section, we will review some known theorems that may be used to 
establish square-function estimates on singular integral operators. 

We begin with one of the hrst such results, the Christ-Journe Tl theorem from 
Section 2 of [C.I87| , which is a square function analogue of the well-known result of 
David and Journe [D,I84| . 
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Theorem 5.1. Suppose that the family of linear operators {0t}t>o given hy 

<dtf{x)= [ tfjt{x,y) f{y)dy 
for some kernels tpt that satisfy 

\Mx,y)\ < Co 


(<+ \x-y\)^+^’ 

c\y- 


\ll)t{x,y) -'ll)t[x,z)\ < Cl \y- A< \{t+ \x-y\) 

for some constants Co, Ci and some e > 0. If 


I ri{Q) 

\Q\Jo 


QCR"- \Q\ Jo Jq 


|0a(x)i 


2 dx dt 


< C 2 


then we have the bound 


Jo IR" ^ 


where C depends only on the constants e, Co, Ci, C 2 and the dimension n + 1. 


We will use this theorem directly in Section [10.21 below. However, this theorem 
is too restrictive to apply to the operators ©f and ©f. (In particular, Of and ©f 
lack smooth kernels.) There are many generalizations of this theorem; we will need 
the following T1 and Tb theorems from IGdlHHI . (We will define a CLP family in 
Section B) 

Theorem 5.2 f [GdlHHl Theorem 4.5]). Consider a family of operators {0t}t>o 
taking values in P > 0, so that Ot = (Of Of , 0(~'"^), where each 0{ acts 

on scalar-valued L^(]R"), and where for g = (^ 1 ,^ 2 , ■ • ■ ;5p+i) G 1 —>■ 

we set 

p+i 

^tg^'^Oigj. 

f=i 

Suppose that there is some 0 > 0 and some C > 0 such that, for all dyadic 
cubes Q, all integers j > 0, and all functions gj € L^(Aj(Q)), where Aj(Q) is as 
in formula (1^ . we have the estimate 

(5.3) \\Ot(lA,iQ)9f\\L^iQ) < ^f^(Q) < t < 2I(Q). 

Suppose further that for some 0 > 0, some CLP family of operators Qs, and some 
subspace H o/L^(]R"), we have that 

(5.4) ||0tQs^||L2(R„) < ||^||L 2 (Rn) for all h € H and all s <t. 

Finally, suppose that 

r^(Q) r clx dt 

(5.5) / \Otl{x)f^<Co\Q\ 

Jo Jq ‘ 

where 1 denotes the (p+ 1) x (p+ 1) identity matrix. Equivalently, we may require 
that 
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for each 1 < j < p + 1, where 1 denotes the function that is one everywhere. 

Then for all f G H, we have that 

(5.6) r [ |0*/(ai)p^<C||/||i.(K„). 

Jo J ^ 

Remark 5.7. The uniform bound 

(5.8) sup||0tg||i2(R„) < C||5||L2(Rn) 

t>o 

follows from the bound (lOl) by summing over dyadic cubes Q of side-length 2^ ^ 
2^ < t < 2^'^^. In particular, establishing the bound (15.31) suffices to show that 0t 
is defined on L^(K”). 

The major advantage of Theorem [5^ over Theorem l5.ll from our perspective, is 
that we need not have pointwise estimates on the kernels of our operators 0t. Rough 
kernels appear in the theory of second-order equations (see [GdlHHl Section 3]) and 
are an essential part of our treatment of layer potentials for higher-order equations. 
On the other hand, we note that the proof of Theorem 15.21 is an easy modification 
of that of Theorem 15.11 to the best of our knowledge, Gran de la Herran and 
Hofmann are the first authors to treat square-function estimates via Tb theorems 
in this generality. 

We now outline the bounds that we will prove using Theorem 15.21 
In Section ini we will show that and 0f satisfy the estimate (15.3p . Notice 
that of is required to satisfy this estimate for all g € L^(R"'), not only all g G 
WA^ |(R") C I—>■ C^); to deal with this technical requirement, in Section I5TT] 

we will extend Of to an operator defined on all of L^(]R"). This extension is fairly 
artificial and is used only for this technical requirement; a similar extension will be 
used for another purpose in Section [10.21 

In Section[7]we will show that Of satishes the estimate (15.41) for all h G L^(]R"), 
and that Of satishes the estimate (15.4p for all h G |(R”). (We do not need 

to extend this estimate to all h G L^(R").) 

Finally, in Section [U we will show that if 

(5.9) Off(x) = Of(fe±)(x), 

(5.10) Off(x) = ^ Of(fje^)(x), 

7„+i<|7|=m-l 

then of i = 0 for almost every x and t. Thus, the estimate dSH) is valid for 
Ot = of . Indeed, one can see that the splitting (I5.9I) - (I5.10I) corresponds to the 
case when 7 = (0,..., m — 1) in formula (14.21) (that is, all derivatives under the 
integral are in t, s) and the case when each term of the integrand has at least one 
y-derivative, respectively. 

We will not be able to show directly that 0^1 or 0^*1 satisfy the bound (|5.5I) . 
In Section [HI we will show that if 2m > n, then 

0^1(0;) = T 4 (a;) -I- Ofa{x) 

where Tj is a Garleson measure (that is, satisfies the estimate (15.51) 1 and where a{x) 
is a uniformly bounded function. Standard techniques will allow us to control Ofd 
by 0fi, using only the fact that d is bounded (that is, without using any special 
cancellation properties); see Lemmabelow. Thus, a bound on 0^1 together 
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with the equation 0f 1 = 0 will give us a bound on 0^*1 and thus allow us to use 
Theorem 15.21 

However, this argument does require control on 0^, and Theorem 15.21 will not 
suffice to bound 0^. 

We will bound 0^ (giving us a bound on Of) using the following theorem, with 
0r' = 0^ and 0( = (0f,0f). 


Theorem 5.11 f |GdlHHl Theorem 2.13]). Consider a family {0t}t>o of operators 
taking values in 0^+^, p > 0, so that 0* = (0(,0j^^) = (0^ , 0^ ,..., 0f~*’^), 
where each Qj acts on scalar-valued and where for q = {q'qn+i) = 

( 51 , 52 , ■ • -, 5 ^+ 1 ) G T^(R" H> C^+i), we set 

p+i p 

Qtg = ^&lgj, Q'^g'= ^Qlgj. 

1=1 1=1 


Suppose that 0* satisfies the bound dOl), that 0f^^ satisfies the bound ([531) for 
all h € and that there is some subspace H' C 1 —1 C^) such that 0( 

satisfies the bound (15.41) for all h £ H'. 

We define the C, 6-norm as 


Q\ is Jc 


ITIIc .5 = sup 

e{Q)>5 \Q\ 


IT(a;)( 


, dx dt 


Suppose that 

(5.12) ||0tl||c,(5 < Gi + C'i||0f''’^l||c,5 for all S > 0 small enough 


where 1 denotes either the p x p identity matrix or the number one. 

Suppose that for each dyadic cube Q C M", we have a measure /rg such that 

(5.13) dp,Q=(j)Qdx, ||V((>Q||icx,(jin) < Co^((5)"\ onQ. 

Suppose further that for each such Q there exists a vector-valued function bq = 
{b'Q,bP^^) G X H' such that 


(5.14) 

r^(Q) f . drdf 

/ / |0*bg(x)P^<Co|Q|, 

Jo Jq f 


(5.15) 

Jr'^ 

bq{x)\Ux < Co\Q\, 


(5.16) 

Rej 

f dgq > cr, 

Q 


(5.17) 


f b'q dgq < pa, 

JQ 

P < l/(2Ci + 4) 


Then for all f £ H' x L^(]R"), 

(5.18) r [ \Qtf{x)f^ < 

Jo 5r" 1 


This theorem is a local Tb theorem; that is, we may test Qtbq near Q, for some 
bq adapted to our particular cube Q, rather than testing 0(1 in an arbitrary cube. 
There is an extensive body of work devoted to generalizing T1 theorems to Tb 
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theorems and local Tb theorems; see, for example, the survey paper [HoflO] , and in 
particular [MM851 ID,TS85l ISemQOl IChr90) for a few of the important milestones of 
the theory. 

As mentioned above, we will establish the bound (15.121) in Section [HI for = 

0^ and 0( = (0f’,0f ), provided 2m > n. We will construct the measure /rg 
and test functions bg = {bg^bg) at the beginning of Section [TOl and therein will 
establish the estimates (15.14^ : we will establish the bounds (j5.15|) . (I5.16p and (15.1711 
in Sections 110.11 and 110.21 The assumption 2m > n will be useful in Section nni as 
well as Sectionini This will allow us to bound 0^, and so together with Lemma lT5l 
will complete the proof of Theorem 11.61 in the case 2m > n. We will extend to the 
case 2m < n in Section El 


6. The decay estimate (15.31) 


In this section, we will show that the operators 0f and satisfy the bound (15.31) 
for all g-l in L'^(Aj(Q)). 

By formula (14.21) for Of, 



/ dp~^''dlsE^(x,t,y,0)gy(y)dy 


2 

dx. 


By Holder’s inequality 

[ \Qf9^^<C\\g\\LHMQ)) [ [ \dr+'^y'f^-^E^{x,t,y,0)\^dydx. 

JO JO JA^(O) 


Finally, by the bound (j3.3l) on the fundamental solution, 

J Q 

Thus, if k is large enough then the operator Of satisfies bound (15.3p . 


Remark 6.1. Suppose that {z) = dzihj{z) for some 1 < z < n and some 
supported in Aj^i{Q). Then 



E / dr+>^dl,E\x,t,y,0)dy,h^:^iy)dy 


2 

dx. 


Integrating by parts in yi, and applying the bound (13.3p . we see that 

In particular, for any h G we have the uniform estimate 

(6.2) ||0f(9*/l)||L2(Rn) < y||/!,||i2(R„), I < Z < 71. 

This formula will be useful in Section [T] 


We now wish to show that Of satisfies the decay estimate (lOI) . that is, that 
||0f/^'|| l=(q) < for all e(Q) < t < 2£(Q) 

for all fJ G L^(]R”) and supported in Aj(Q). 
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Choose some dyadic cube Q and some t with £{Q) <t< 21{Q). Suppose first 
that S 2), where 2) is as in Remark [2.371 and is supported in Aj^i{Q). Recall 
that (and thus is defined in terms of extensions of fF Thus, 

we begin by choosing an appropriate extension. Let be the function given by 
Lemma [331 we then have that 

Tr- = /^ sup||V'"iJ^'(-,i)l|L^(R") <C||/'I|l=(R"). 

” t<o 

We may assume without loss of generality that TrF-^ = 0 outside of 2-^+^Q. Let 
r]j{x,t) be smooth and satisfy the bound |V*77j(a;, t)| < Ci{ 2 H{Q))~^, with 

rij{x,t) = 1 if a; £ 2^'^'^Q and — 2H[Q) <t < 2H{Q), 
rjj{x,t) = 0 if a; ^ 2 ^^^Q or |t| > 2 ^^^^{Q). 


Then | F[^ = fF We take = rjj HF Observe that if j > 2, 

then V™F^(a;,<) = V^H^{x,t) = 0 if |t| < dist(a;,]R" \ 2^-^Q). Furthermore, we 
still have the bound 




Now, by the definition (14.31) of Of and by formulas (12.3211 . and (12.3111 




2 k 


E 

\P\=m 


dr^^dl,E^{x, t, y, s) {AV^F^ ) 0 iy, s) ds dy 


dx. 


Applying the bound (|3.3I1 . we see that 

<C2-^-(2fe+0||/a ||2^^^^) 

for all G 2) supported in Aj^i{Q); by density we may extend to all f G | (R") 

supported in Aj^i{Q). 


6.1. Extending Of to all of L^(K"). We now must extend Of to an operator 
defined on all of that still satisfies the estimate (lOll . Essentially, this 

argument consists of defining a projection operator from to WA^ |([^")) 

the space on which Of naturally acts. 

Because L^(]R”) is a Hilbert space, there is an orthogonal projection operator 
Ow ■ !->■ 1FA^|(R"'). For example, if m = 1 then Owf = 27||W, where 

A||U = V|| ■ f. This is the most natural mapping from L^(]R”) to |(R"); 

however, this mapping does not satisfy adequate decay estimates. Thus, we must 
refine this mapping by applying cutoffs before and after projecting. 

Let Wj be the closure in L^(]R”) of 

{l2i+2Q TCm-l + (1 ~ 12J+2 q)/ : ‘P G (7“, / G L^(]R"')}. 

Loosely, elements of Wj are higher-order traces in the cube 2 W‘^Q and are merely 
arbitrary arrays outside of that cube. Let Oj denote orthogonal projection 
from L^(R”) onto the subspace Wj; observe that Ojf = f outside of 2W'^Q, 
Furthermore, if (/? is a nice function then OjifTrm-i f) = Tcm-i P- 
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Let r]j be a smooth partition of unity; that is, = 1 for t near zero, 

with rjj supported in x and satisfying < 

C2-^HiQ)-^ for all {x,t) e 

Define ttj : Wj i-A as follows. Suppose that / = Trm,| f in 

for some smooth function (f. We may renormalize (p so that Jg Tr d''(p = 0 for all 
1^1 < m — 1. Let TTjf = Tr^ |(? 7 j(^). We remark that iTjf is well-defined, that is, 
^m,\(,r]jV’) depends only on Trm,| ‘P- Furthermore, observe that TTjf is supported 
in Aj^i{Q). Finally, by the Poincare inequality 

lkj/||L2(A,-,i(Q)) < C'2'^"^^||Tl'm,| <P||l 2(2I+2Q) = (72-^”/^ || /1| L2 (2i+2Q). 

We will extend to an operator on all of using the orthogonal projection 

operators Oj. Observe, first, that Tt'jOjf = 0 outside of Aj^i{Q), and second, that 
if / = 0 in then Ojf = f and so TTjOjf = 0. 

If ip is smooth and compactly supported, and renormalized as above, then 

OO OO OO 

TFm.l P = ^ TTj (Oj (TF„_i (/?)). 

j=0 j=0 j=0 

We define Of*/ = f)- Now, if is supported in Aj(Q), ob¬ 

serve that TTiOif = 0 for all* < j — 2, and furthermore that ||7riOi/'Jj|i2(Rn) < 
C'2™/2||/J|li2(^^(g)) for all i > j - 1. Then 

OO OO 

i=j-l i=j-l 

OO 

i=j-l 

and so Of* satisfies the decay estimate (15.31) provided k is large enough. 


7 . The quasi-orthogonality estimate (15.41) 

A family of operators {Qs}s>o is defined to be a Calderon-Littlewood-Paley 
family, or CLP family, if 

Qsf{x) = [ s"” (p{y/ s)f{x- y) dy 

JR" 

for some (p G L^{MA) that satisfies the conditions 

(7.1) 1^(01 <Cmin(|C|Me|--), \p>{x)\<C{l + \x\)—- 

for some cr > 0, where ip denotes the Fourier transform of p, and such that Q 
satisfies the conditions 

(7.2) \\Qsf\\L^(R^) + ||sV(5s/||L2(Rn) < C'||/||i 2 (Rn) for all s > 0, 

(7.3) [ riQ./(cr)P^<C'||/||i.(«„), 

J M”' « 0 ^ 

(7.4) rQi-=i 

Jo S 

where convergence to the identity in the last formula is in the strong operator 
topology on S(L^(IR")). 
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We remark that, up to multiplying by a constant, it suffices to require that 
kp satisfy the bounds FTTl the bound |^(^)| < be radial, and be real-valued 

and not identically zero. 

Specifically, notice that 

&o = ^isom- 


Given this relation, the estimates (IL2I) and dm) follow from the estimate |<^('C)| < 
Cmind^l*^, 1^1“^) by Plancherel’s theorem. To establish the identity (17.41) . we nor¬ 
malize ip as follows. Because ip is real-valued, so is p. Thus, the integral 

nOC 1 

/ p{sO^-ds 
Jo ^ 

is independent of ^ provided ^ ^ 0. If (/? is both radial and real-valued, then p 
is real-valued, and so this integral is positive. We may normalize p so that this 
integral equals 1. It is then straightforward to establish the condition dH. 

In this section, let p be bounded, radial and supported in B{0,2) \ B{0, 1/2). 
(In Section 19.21 we will use a CLP family again; in that section it will be more 
convenient to take p, rather than p^ compactly supported.) We wish to establish 
the bound (|5.4|) . for the operators 0t = 0f or 0* = 0/^. We proceed as in [GdlHHj . 

We begin with 0f. Fix some h € L^(R"'); we wish to bound QfQsh. For each 
1 < j < ’T') let // satisfy 


m) = 




27rz|^|- 


:Qshj{0 = 


2Tri\^\‘‘ 


■p{sO /17(C)- 


Then Qgh^ = and so 


0fQ,Mx) = ^0f(a,/^)(x). 

t=i 


By the bound (j6.2|) . 


n 1 ^ 

lefQshWmm^) < ^||0f(a,/J')llL=(RG < 




t=i 


Notice that |^(C)| < C|C|, and so ||/||L 2 (Rn) < Cs\\h\\L 2 (^R„y Thus, 

||0f(5s^||L2(R") < G-^||/l||L2(Rn). 

Therefore, 0f satisfies the bound (15.41) for 0 = 1 (and thus for any 9 < 1). 

We now consider Qf‘. Let the subspace B be WA^ |(R"'); recall that this is the 
natural space upon which and 0f act. It suffices to show that 

|| 0 fQs Tr,„,| pWl^r^) < ‘P\\Li^(R^) 

for some 9 > 0 and for all smooth, compactly supported functions p. To establish 
this bound, we begin with the following lemma. 


Lemma 7.5. If p is smooth and compactly supported, then 

Qs{ThCm,\ p) = Tr^.l 
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for some function that satisfies 

Proof. For each j with 0 < j < to — 1, let f^{x) = 0). Observe that if f3 

is a multiindex and /3„+i < \j3\ = to, then 

p)p = dx" 

Qs is a convolution operator and so commutes with horizontal derivatives, and thus 
(Qs Trm,| = dx"Qsf^’'+^{x,0). 

Let gi = QsP. Then gl{f) = (p{sf) P{f) and so < C's|/^(^)|2. Thus, 

[ \9iiOW^-^^-^d^<Cs [ \P{OW^-^^d^ = Cs [ \V^-^P{x)\^dx. 
Jr-^ Jr^ Jr^ 

Because P\ < |Trm,| we have that 

[ df < CsIlTh^,! 

Jr" 

Let Qs = (5s(Trm-i <f)', then satisfies (gs)^ = for each I 7 I = to— 1. We 

have established that 

JR’' 

Extending Qs using Lemma 13.41 completes the proof. Notice that more general 
extension theorems, also appropriate in this case, are well known; see, for example, 
[Tri78l Theorem 2.9.3] or |Tri83[ Theorem 2.7.2]. □ 

The estimate (15.41) for Of follows quickly from Lemma 17.51 By the defini¬ 
tion (14.3p . 

ef(Q, i^m,\ P){X) = ^s){x, t). 

But by the definition (12.301) . 

0 f(Q,T>^,|(p)(x) =-Car+"n^(l-dlV™4>,)(a;,t). 

But u = n^(l_ 74 V'"<i)s) satisfies Lu = 0 in furthermore, by the bound (12.151) . 

IV'"'“IIl 2 (r"+ 1 ) < C'IV™$s|li 2 (Kn+i) < CyillTF^j (^||i 2 (Rn). Applying the Cac- 
cioppoli inequality and Lemma 13.21 in small cubes of sidelength t/C suffices to 
establish that Of satisfies the bound (15.41) for 0 = 1/2. 

8. The semi-horizontal single layer potential 

In this section we will prove the following theorem. 

Theorem 8.1. Suppose that^n+i < ItI = to— 1. Then we have the square-function 
estimate 

[ \Qf{geP{x)f^dxdt < C\\g\\l2^R,^y 
7r"+i I 










32 


ARIEL BARTON, STEVE HOFMANN, AND SVITLANA MAYBORODA 


Proof. Let Qtg = 0f (ge^). We want to apply Theorem l5.21 As shown in Sections[ 6 ] 
and [7l the bounds (15.31) and (15.dp are valid for this choice of 0t. We are left with 
the estimate (|5.5I) . 

Recall that by formula (14.21) . 

Otgix) =t^ [ t, y, 0) g{y) dy. 

In particular, 

&tl{x)=t>^ [ dr+’^dl,E^{x,t,y,0)dy. 

Js." 

Let j satisfy 1 < j < n and 7 j > 0; by assumption on 7 such a j exists. Let 
C = J — Bj + Bn+i- By formula (I2.26p . we have that 

0 a(x) = [ dy^{dr+’^-^dl,E^ix,t,y,o)) dy. 

By the bound (13.3p . for almost every {x,t) G if k is large enough then 

v{y) = df^+'^-^dl,E^{x,t,y,0) lies in both L\M.^) and in Thus, 

[ dy, [dr+’^-^d^^E^x, t, y, 0)) dy = 0 

JR^ 

for almost every {x,t) G Thus, 0*1 = 0, and in particular the bound (15.51) is 

valid. □ 

We are now left with the double layer potential Of and the vertical single layer 
potential 0^/ = Qf{fe±). In the following sections we will use the full force of 
Theorem 15.111 to bound 0^; we will bound as a component of the auxiliary 
operator 0^ in Theorem 15.Ill 

9. The Carleson estimate (I5.12p 

We will let 0t = (0^*, 0f), with 0^ denoting the purely vertical component of 
0f (that is, = 0f(/e_L)). The bounds (15.31) and (15.41) are valid. We wish to 
show that the bound (15.121) is also valid. Recall that this bound is given by 

|| 0 tl||c ,<5 <Ci+ C'i|| 0 ^ 1 ||c ,5 if d > 0 small enough. 

In Section m we showed that 0f(le.y) = 0 whenever 7 „+i < I 7 I = m — 1; the 
bound (I5.12p . with 0( = 0f , follows immediately. Thus, we need only bound 0^*1 
by 0^1 and a constant. 

Recall that by formulas (14.31) . (12.321) and (I2.31L if / = Tr„j | F, then 

0f/(a;) = - Y C /" df^+'^d“^E^{x,t,y,.s) Aai 3 iy,s)d^F{y,s)dsdy. 

|a| = |/3|=m 

Using the bound (13.31) . we see that if is bounded then the integral converges 
absolutely for almost every {x,t) G R"^^. 

Recall that 0^* acts on arrays of functions of the form tp = Tr„j | (/?; these arrays 
if are indexed by multiindices /3 with /3n+i < |/3| = "m. Fix some such /3. By 
choosing 
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we see that 
(9.1) OF^isix) 



t, y, s) Aayiy) ds dy. 


We use formula (12.261) to convert one of the derivatives in t into a derivative in s; 
we evaluate the integral ds to see that 



dr+^-^d^^,E^{x,t,y,0)A^p{y)dy. 


Observe that this is a sum of terms depending on a and jd. In Section [^?T] we will 
bound the terms for which a„+i > 0. In Sections [9.2H9.51 we will bound the terms 
for which a„+i = /3n+i = 0 ; this case is the most involved, and will closely parallel 
the argument in |GdlHHI Section 3.1]. Finally, in Section l9^ we will bound the 
remaining terms, that is, the terms for which an+i = 0 and /3n+i > 0 ; this bound 
will rely on the bound in the case a„+i = /3n+i = 0 . 


9.1. Terms with an+i > 0. Observe that if a„+i > 0, then a = 7 + e„+i for a 
unique 7 with I 7 I = m—1. For any such 7 , let 7 = 7 +e„+i. Then by formula (I2.26L 

^ d^^A"'+'^-'E^(x,t,y,0)A^/s(y) = - E dAdr^^d^^(x,t,y,0)A^/s(y). 

|Q;|=m | 7 |=m—1 

ck 7 i+i >0 

Thus, 


of e^(x) 


[ dr+'^-^dAE^{x,t,y,0)A^p{y)dy 


- dA^dl,E^{x,t,y,0)A:^is{y)dy. 

Jr- 


In this section we will bound the second sum; we will consider the first sum in 
Sections 19.2119.61 

By formula (14.21 . the second sum is equal to Of a/ 3 , where (ap)^ = Aj^p. Notice 
that ap is bounded. We may control Qfap using a standard technique in the study 
of T1 theorems. Let 

Ptf{x) = f(y) 

where "0 is smooth, nonnegative and satisfies fg„ 0 = 1. We do not require that 0 
be compactly supported. Fix some 7 with I 7 I = m — 1, and let 

'l'ta(a;) = Of {aA){x) — Pta{x) Qfe^{x). 


Then 

||Pt(a/3).y0fe^||c,5 < ||a||L“||0fe^||c,5. 

Either 7 = 7 _l and so Ofe..^ = 0f 1, or 7^+1 < I 7 I and so Ofe..^ = 0. (See the 
proof of Theorem 18.11 1 So to bound || 0 fd/ 3 ||c, 5 , we need only control 'I'(a(ai) for 
arbitrary bounded functions a. 

We wish to show that Ijiftallc = ||4ita||c,o < C'llajlioo. (We will not need the 
|| 0 fl||c .5 term on the right-hand side; thus, we will simply bound the left-hand 
side for J = 0 and observe that this gives our desired bound for any positive 5.) 
We will do this by applying Theorem 15.21 to Observe that 4itl(x) = 0, and so 
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the bound (15.5p is valid. We need only verify the bounds (|5.3I1 and (j5.4l) for the 
operator dtj; because these bounds have been verified for 0 f, we need only consider 
Tta(a;) = Pta{x) Qfej{x). 

We begin with the bound (15.41) . Observe that 

i^h{o = ${to 

where Qs is the operator defined in Section[71 Recall that ^ is supported in B{0, 2)\ 
B{0, 1/2). If we require that be smooth and supported in B{0, 1/2), then PtQgh = 
0 whenever s < t; thus Tj satisfies the bound (15.41) . 

We now establish the bound (15.31) . By Section [51 we know that the operator 0f 
satisfies the decay estimate dOl). From this we may verify that, if Q C K” is a 
cube with £(Q) <t< 2£{Q), then 

If Tp is smooth as well as being supported in B(0, 1/2), then ip is a Schwartz function 
and satisfies the estimate \ip{y)\ < C'Ar(l + |?/|)“^ for any N > Q. If gj is supported 
in Aj{Q), and £{Q) <t< 2£{Q), then 

sup|P*g,(2^)l < CNt-2-^^\\g,hi^AAQ)) < 

xeQ 

and so if we choose N large enough, then the operator satisfies the bound (j5.3p . 

Thus, by Theorem l5.2l we have that the operators ikt satisfy the square-function 
estimate (lOl) . To show thcit ||TiCz||(;; ^ ^11 ^11 (R ”^)7 need only show that tho 

estimate dSH) for test functions implies an estimate for L°° test functions. 

Lemma 9.2. Suppose that the operators satisfy the square-function estimate 

r [ |vI/,/(cr)| 2 ^<Co||/||i.(r) 

Jo iR" ^ 

for all f G i—> C), and that for some 9 > —2, dtj satisfies the off-diagonal 

decay estimate 

\\'^tg,\\L^^Q) < Ci2-^-("+2+«)/2||5,|U. if£(Q) < t < 2£{Q) 

for all j > 1 and all gj supported in Aj{Q) = 2l+^Q \ 2‘^Q- 

Then there is some C depending only on Cq, Ci and 0 such that ikt satisfies the 
Carleson condition 

1 f ,, \\ 2 dtdx 2 

TmJjo s c'RiL-c-) 

for all bounded functions b. 

Proof. Choose some cube Q and some bounded function b. Let bj = 61^ . (g^; recall 
that bo = ^Iao(Q) = bl 2 Q. Then by the square-function estimate, 

/ r r^iQ) 

[JqJo <qi&0||L=(2Q) <C||6|U^(Kn)|Qr/^ 

Furthermore, if j > 1, then by the decay estimate applied in cubes R G Q with 
side-length t < £{R) < 2t, 

/ rRQ) r m \ 
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Summing in j, we see that 
r^(Q) 


Q Jo 


dtdx \ 


iQr/" 


as desired. 


□ 


9.2. Terms with a„+i = fin+i = 0: preliminaries. We have now established 
that 

Qfefi{x)= XI / d'^""'"~^d^,sE^{x,t,y,^)Aafi{y)dy-Qfap{x) 

. _r. 


Q!n + 1 —0 


and that 0f d /3 satisfies the desired estimate. In order to show that the remaining 
sum satisfies the estimate (15.121) . we need only show that the operator ©f, defined 
as 

(9.3) 0f/(x)= X / 5™+''"^a“i;^(a;,t,2/,0)Ta/3(2/)/(2/)d2/, 

a„+i =0 

satisfies the Carleson-measure estimate 

(9.4) ||0fl||c.5<Ci+Ci||0,^l||c,5 

for any multiindex /3 with /3n+i < |/3| = m. 

We will make use of the horizontal operator L||, defined as follows. Recall that 
L is an operator acting on /Q^(M"+^)-functions. We may (formally) define the 
operator L||, acting on ;Q^(K”)-functions, as 


(9.5) 


4^11/= (-l)^ 




|a| = |/3|=m 

;n + l=/3n+l=0 


where 5“, are understood to be derivatives in the n horizontal directions. (The 
operator L|| has a weak formulation, as in formula (1131).) 

To establish the bound (19.41) for (dn+i = 0, we will follow the argument of |GdlHHl 
Section 3.1]. We remark that the argument we will make in this section is valid only 
in the case where the order 2m of L (and thus Ly) satisfies the inequality 2m > n. 
Thus, the argument of Sections [9] and [TO] will only establish boundedness of Of 
and 0^ in the case of operators of very high order. In Section [11] we will show 
that bounds on Of and Of, for operators of high order, imply the corresponding 
bounds for operators of lower order, completing the proof of Theorem 11.61 

We will use some tools from the proof of the Kato conjecture, in particular from 
the paper [AHMTOT] . The following lemma was established therein. 

Lemma 9.6. Suppose that 2m > n. There is some W depending only on the 
standard constants such that, for each cube Q C M", there exist W functions fQ,w 
that satisfy the estimates 

(9.7) I |V]p fQ,wf < C\Q\ for any cube R with £{R) = £{Q), 

c 


|i||/Q,^(ai)| < 


(9.8) 
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and such that, for any array 74 , 


(9.9) 


I r^Q) 

TW\Jo 



dx dt 
t 


< 


/ / \('ytix),AfVffQ^Ux))\ 

Q \Q\ Jo Jq 


, dx dt 


where Aff{x) = fg, f{y)dy, for Q' C Q the unique dyadic subcube that satisfies 
X G Q' and t < i{Q') < 2t. 

Specifically, the bound (19.81) is the bound (2.19) in [AHMTOl] . The bound (19.71) 
follows from the bound (2.18) in |AHMT0T| (if R = Q) and the observation that, 
by Lemma 3.1 in [AHMTOl] and the definition of fq^w therein, fq^w = ^'^fR,w 
whenever 1{Q) = t{R). Finally, the bound (19.91) is simply Lemma 2.2 of [AHMTOT] . 
The requirement that 2m > n is a sufficient condition (see [AHMTOTI Proposi- 
ton 2.5] or [Dav95[lAT98| l for L|| to satisfy a pointwise upper bound; this condition 
is assumed in the proofs of the above results. 

Let ( 7 t )/3 = l 5 <t<i/ 50 f 1; notice that 74 includes ©f 1 provided Pn+i = 0. Thus, 
the estimates (19.41) and thus (15.121) . for fJn+i = 0, follow from bounds on the 
quantity 


1 


I 


min(l/( 5 ,^(Q)) 


I ^tHx)Afd^^fq,^,{x) 

/3„+i=0 


“^dxdt 

t 


for 6 < £{Q) and d < 1. 

We will divide this quantity into a sum of controllable terms as follows. Let 
Ptf{x) = f * iftix), where = t~'^ij;{x/t). We require that ijj be smooth and 

nonnegative, that J ip = I, and that il^{x) = 0 whenever \x\ > 1/2. (We will later 
impose some additional constraints on i/;. Notice that it is convenient to use a 
different approximate identity P 4 in this section from that used in Section [9.Il l Let 

R\'^F{x) = ©f l(a;) (A?F(x) - PtF{x)), 

R^/Pix) = e^lix)PtFix) - e^iPtF)ix), 

R^/Fix) = e^iPtF-F)ix), 

Rf^F{x) = ©f F(cr) 


so that we seek to bound 


/3„+i=0 i=l /3„+i=0 

We begin with Observe that by the definition (19.31) of ©f, 


^ = X ^E^{x,t,y,0)Aap{y)d^fq^^u{y)dy. 

.+ 1=0 a„+i=0 

hn + l=0 
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Using formula (12.2611 and then integrating by parts in y, we see that this quantity 
is equal to 


/3„+i=0 


f drdt^E^{x,t,y,0)^fQ^Uy)dy 


and by the bounds (19.811 and (13.311 . we see that if k is large enough then 



‘^dxdt 

t 


< C. 


9.3. The term R\’^d^fg^w Next, we bound R\’^d^fq^w We begin by bounding 
©f l(a:). Recall the following special case of Morrey’s inequality (see, for example, 
[Eva981 Section 5.6.3]): if a: G Q C K", then 

m , . 1/2 

|T(a;)| < ^ C£{Qy ( -r |V||Up j provided 2m > n. 
i=o ' 

We apply this bound to the function v{x) = (x,t,y,0), a locally 

Sobolev function for almost all y and t. Then by the bound (13.3L we have that if 
|a| = m and either \t\ = £{R) or |t| < 1{R) and j > 1, then 


(9.10) [ |a“,9r+'=-^U^(a:,<,j/,0)|2dy < C^(i?)-”-^'=2-^(”+2'=). 

Jaj{r) 

Observe that this bound is valid for any k > 1I2 — n 12] in this section we will need 
this bound only for k large, but in Section [TUI we will need this bound for fc = 0 and 
fc = 1 as well. Also, by formula (12.211) it is valid with the roles of y and x reversed. 
Using Holder’s inequality and summing over j, we see that if A; > 0 then 


|0fl(x)| 


^ ^E^{x,t,y,0)Aai3{y)dy 


< C 


and so 

dl^fQ,w{x)\ < ClAf ( 9 %,„(x) - Ptd^fQ^^{x)\. 

Thus, we need only bound Af d^fq^wix) — Ptd^fq^w{x). We will do this using a 
standard orthogonality argument. 

Let R\ = Af — Pt- Recall that the kernel ip of Pt is supported in 5(0,1 /2); thus, 
if X G Q and t < £{Q), then R]F{x) = (12 qU)(x), and so 

f \RyPd%,^ix)\^^ r f (l 2 Qa%,.)(x)p^. 

Jq >' Jo Jr^ ^ 

Let {Qs} be a CLP family, as in Section 0 but with the kernel ip (and not its 
Fourier transform p) supported in 5(0,1/2). By the identity (17.41) . 



rRQ) 

Jo 


f \RyPd%,^{x)y^<c 

Q JO 


nOO p pOO 

/ / / -RtQ?(l2Q9%,„)(x) 

Jo 7k" Jo 


ds 


■ dx dt 
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By Holder’s inequality, for any number £ > 0, 

We claim that 

/s 

(9.11) ||i?t^Qs 5 '||L 2 (K.) < Cminf|| 5 ||i 2 (K„). 

Choose e = 1/6. Assuming validity of the bound (I9.11|l . we have that 


Interchanging the order of integration and evaluating the integral in t, we see that 

Jo JQ ^ Jo iR" S 

By the bounds (17.311 and (j9.7ll . we have that 


pf-iQ) r „ iJt df 

/ / < c||a%,„||i2(2Q) < c|Q| 

^0 Jq ^ 

as desired. Thus, to complete our bound on Rl’^d^fg^w, we need only establish the 
estimate (19.111) . 

Suppose first that t < s and so mm{s/t,t/s) = t/s. By definition of Rf and Qs, 


1 q'(2 /) -tptix- y)j ifsiy - z) g{z) dy dz. 

Notice that y^lg/ (y) — iptix — y) dy = 0, and that the integrand is zero unless 
\x — y\ < Ct. Thus 

R\Qsg{x)< j I (-^lQ,{y)-'ijjt{x-y)]{ips{y-z)-(ps{x-z))dyg{z)dz. 

JR’^ JB{x,Ct) \N \ J 

Suppose that y G B{x,Ct). Because ip is supported in 5(0,1/2) and s > t, if 
|a; — z| > 2Cs, then ipsijj — z) — ips{x — z) = 0. Otherwise, 

{‘Psiy - z)- ips{x - z)\ < Cs~'^~^\y - x\. 




RlQsg{x) = I / ( 


Thus, 


RlQsgix) = C 


Ib(x,2Cs) Ie 


B(x,Ct) 


—^lQ,{y) - ipt^x - y) 


dy \g{z)\dz 


<C-i 
s Je 


t 


\g{z)\dz < C-Mg{x) 

B{x,Cs) S 


where Mg denotes the Hardy-Littlewood maximal function oi g. It is well known 
that M is bounded LP(W^) i—>■ LP{W^) for any 1 < p < oo, and so the estimate (19.111) 
is valid whenever t < s. 
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Recall that the kernel (pa of Qa also integrates to zero and that the kernel V't 
of Pt is also smooth. Thus, by a similar argument, if s < t then \\PtQag\\L^{Rn) < 
C'(s/^)ll5llL2(Rn). Bounding pffQsff is somewhat more involved, because the kernel 
]^1 q' of is not smooth. 

Suppose s < t. Let r] = r]t^a,x be a smooth cutoff function that is identically 1 in 
Q' and is supported in (1 + \fsJt)Q' . Then |V? 7 | < Let 

B^G^x) = y riiy) G{y) dy. 

By the same argument as above, we may show that if s < t then \B^^Qag{x)\ < 
C{s/tyGMg(x). To conclude the argument, notice that 

\AfQag{x) - Bf^^^Qag{x)\ <^[ \Qsg\- 

^ Jsuppr}\Q' 

Notice that jsupp?? \ Q'\ < CV^\Jsjt. We apply Holder’s inequality to see that 

\A^Qag{x) - B%Q,g{x)\ < Ci^- J j 

/ \ 

<g(^jJ M(|Q,g|3/2)(a;)2/3. 

Because 3/2 < 2, the estimate (19.lip is valid for s < t as well as t < s. This 
establishes our desired bound on d^fg^w 

9.4. The ter m Next, we bound We will use the following 

lemma from [AAA~*~ll] : this is a square-function Tl theorem that is somewhat 
simpler than Theorem 15.21 but has more stringent requirements. 

Lemma 9.12 (Lemna 3.5(ii) in |AAA~*~ll] '). Suppose that {Rt}t>o is a family of 
operators satisfying 

(9.13) IIA(ri.4,(Q,)IIL,«, < 

for all 0 < t < £{Q) and all j > 1. Suppose further that for all t > 0, all F G 
L^(]R”), and all well-behaved veetor-valued funetions F, we have the bounds 

(9.14) ||Rt^’||L 2 (R") < C'||^"||l 2 (r"), ||Rt div|| ,F||l 2 (r") < y ||.F||L 2 (Rn). 

Finally, suppose that i?*! = 0 for all t > 0. 

Then 

JR-^ Jo t 

for all F GL^{W^). 

By the definitions of 0f and Rf’^, we have that 
Rf^F{x) = Q^,l{x)PtF{x) - Q^{PtF){x) 

= E / dr+^-^d^E^{x,t,y,0)A^f,{y){PtFix)-PtF{y))dy. 

a„+i=0 
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Let 

(9.15) =t^ [ d^dT+^-^E^{x, t, y, 0) A^p{y) {PtF[x) - PtF{y)) dy. 

Jmv- 

In this section we need only bound for a„+i = 0; in Section 1^75] we will need 

an estimate on in the case where a„+i > 0. 

Observe that Ptl{x) = Ptl{y) = 1 and so — 0. 

Now, recall that PtF{x) = f t~^ip{{x—y)/t) F{y) dy for some smooth, compactly 
supported function ip; then 

\\PtF\\L2(^n^ < C'||F||i2(Rn) and ||Pt(diV|| F)||i2(Rn) < y ||F||i2(Rn). 

We must use this fact to establish the bounds (19.131) and (19.141) . 

Let Q be a cube and let 0 < t < i{Q). If j > 1 and F is supported in Aj{Q), 
observe that PtF{x) is supported in Aj^i{Q). Thus, by the bound (13.311 . 


(In the case j = 1 some extra care must be taken to establish this estimate; however, 
it may be done by considering the cases t > i{Q)l2 and t < separately.) 

This implies the bound (19.131) . We are left with the uniform I? bounds (19.141) . 
Suppose that F is supported in 8Q and that £((5)/2 <i< Then PtF{y) = 
0 for all y ^ 16Q and so 


|i??’“’^n^)| < Ct’^ [ \d!^dr+'^-^E^ixA,y,0)\\PtFix) - PtF{y)\dy 

<ct'^ f - PtFiy)\dy 

Jl6Q 

oo y 

+ Ct'=^|PtF(a;)| / |5“5r+"-'P'^(x,t,2/,0)|dy. 

Jah(Q) 


Applying the bound (19.1011 . we see that 

|P?’“’^F(x)| < C\PtF{x)\ + Ct-"/2||P*P|U2(i6Q). 


Thus, 


||Pt(Pl8Q)||L2(Q) < C'||Pt(Fl8Q)||L2(Rn). 


We sum over cubes of side-length t; this yields the bound 


WRtPhHQ) < C'||PtP|U2(R„) 

and, combined with the existing bounds on PtF and Pt div|| F, yields the desired 
estimates (|9.14ll . 

Thus, Lemma l9. 121 applies and we may bound the operator In particular, 

using the bound (19.711 and arguing as in the proof of Lemma l97^ we have the desired 
Carleson bound on R^’^d^fg^w. 
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9.5. The term R^’^d^fQ^w Finally, we consider the term R^’^d^fg^w As in the 
case of i?tbut unlike R]’^ and , we will not be able to bound the individual 
terms Rf’^d^fg^w', we will only be able to bound 

RffQM^)= E 

/3„+i=0 

a„+i=0 
/3n + l=0 

Another complication of this section is that we will only be able to establish the 
bound 

\\R!fQ,n.\\c,S<C,+Ci\\&^l\\c,S 

and not simply a bound of the form \\RffQ,w\\c < Ci- 

Let / = (/q,uj — PQ,w)vQi where pq^w is an appropriate polynomial of degree 
TO — 1 and where rjq is a smooth cutoff function that is identically 1 in 2Q and is 
supported in iQ. By the bound (|9.7I) . 

||V|p/|U.(M„) < C\Q\^/\ 

Furthermore, f = in 2Q whenever |/3| = to. Using the bound (19.101) on 

and the bound (EH) on V^fq^w, we may show that 

1 r r^(Q) _ _ rlfrlT 

W\ jq I < C 

and so we need only bound R^f{x). 

By the definition (19.51) of L||, and by formula (I2.21L 

R^ix) =t^ [ ar+''-iL|U^*( 2 /, 0 ,x,t) (PJ - f){y) dy 
Jr^ 

where Lj| is taken in the y variable. Recalling that L* {y,s,x,t)) = 0 away 
from (x,t), we see that 

L|U^*(y,0,a:,t) = (-1)™+^ E 0 ))' 

«n + l+Cn + l>l 

Thus, we need only bound the quantities 

R^’^’^fix) = [ diMMdisdr^'^-^E^{^,t,y,0)){PJ - f){y)dy 

Jr'‘ 

where at least one of Cn+i and ^n+i is positive. 

For each multiindex (, we write ( = ^|| +C-Le_L, where C_l = Cn+i and where C|| is 
a multiindex with (C||)„+i = 0. Integrating by parts and applying formula (I2.26L 
we see that 

Jr-^ 

We wish to bound for 0 < t < (.{Q). Let ^ C Q be a dyadic subcube 

with t/2 < ^{S) < t. Let 'dj ^ smooth partition of unity with rjj supported 
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in y4j^i(S') and with the usual bounds on the derivatives of ijj. Let fj = rjj /. Then 

OO « 

j=o 


By the bound (jQ.lOL if cc G 5” then 


' [ t, y, 0) {Ptfj - fj){y) dy 

JR'^ 

< ct-"- 2 a 2 -^("+ 2 '=+ 2 C^) [ laCii (PJ, - /,)(y)|2 dy. 

Ja,,(S) 


Thus, 


-j{n+2k+2C±) 


/A,-2(5) 


p OO 

''Q scQj =0 

Summing carefully, we see that 

[ \R^,’^’^f(x)\^dx<Ct-^<- [ \d<HPtf - f){y)\^ dy. 
JQ 

Now, by PlancherePs theorem, 


\d<nPtfj - fj)iy)\^ dy. 


Jo Jr^ t 

< [ |/(a;)p f {1 - ^ duj 

Jr" Jo t 

where tjjtix) = t~'^'ip{x/t) is the convolution kernel of Pt. We require that ijj be 
radial and make the change of variables s = t|a;|. Then 


\d‘^nPtf-f)iy)\^dy^ 

Jo Jr^ ^ 

< f |..pl^l|7(a2)p r s-^‘^-0-ks)r-du2. 

Jr" Jo '5 

We require that f ip = 1, and that the higher moments are zero, that is, that 
/ x^ ip{x) dx = 0 for all |d| small enough. This implies that 1 — ip{s) is small for s 
small and so — ip{s)Y is integrable near zero. Because ip is smooth and 

compactly supported, we have that ip{s) is bounded. If C-L > 0 then the integral in 
s converges. (We have that ip{s) —>■ 0 as s —>■ oo, and so the integral must diverge 
if Ct = 0.) Because WyffWmm^) < Cy^l we have that 



<C\Q\ 


whenever Ct > 0. 

'~3 c A 

We are left with the terms R^’^f for = 0; recall that we need only consider 
C_L + > 1 and so we may assume Ct > 1- Because C_l = 0, we have that 


Py’V(x) = t’^ [ AMdisdr^’'-^E\x,t,y,0) {d<Ptf{y) - d^f{y))dy. 
Jr^ 
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Recall from Section [93] that the operator given by formula (|9.15|) . satisfies 

square-function estimates. Thus, we need only bound 

= t^ [ A^c(y) t, y, 0) id‘^Ptfix) - a^/(y)) dy. 

Let 7 = ^ — ej_. We use formula fl2.26p ; we then see that 


= [ Aic{y)dlsdr^'^E^{x,t,y,0) {d^Ptf{x) - d^f{y))dy. 

JR" 

We recognize the integrand as being much like the kernel of the single layer poten¬ 
tial. By formula (14.21) for 0f, we have that 

Rf^'^fix) + = -Pt{d^f){x) Ql{A^c^^){x) + 0f d^fe^){x). 

If 7 7 ^ 7 _l, then the operator 0* given by 0t/ = Qf{fe^) satisfies the bound (15.61) 
(see Section IH). If 7 = 7_ l , then 

0f(kI^c5Ve,) = 0,^(kI^c9V). 

The operator 0^ = satisfies the conditions of Theorem l5.2l albeit with 

constants depending on || 0 ^ 1 ||c_ 5 , and so also satisfies the bound (15.6|) . Thus, in 
either case, we have the bound 

r PR rlxdf 

/ / \Qf{A^^d^fe,){x)\^^ < (C + C|| 0 ^||c,^)|Q|. 

JR-^ JS ^ 

To bound Pt{d^ f){x) Of {A^(^e^){x), recall Carleson’s lemma (see, for example, 
[Ste931 Chapter II, Section 2.2]). 


Lemma 9.16. Let P{x,t) be a function and dy be a measure defined on 
Then 


/. 


P{x, t) dy{x, t) 


< c( sup 
VrcR" 


li?l 


U>^{L 


sup \F{y,t)\dx 

\x-y\<t 


provided the right-hand side is finite, where the supremum is taken over cubes R C 
R". 


We wish to bound 

r ^min(l/< 5 ,f(Q)) i 

1 ^ \Pt{d‘^f){x)\^\Qf{AiCe-^)ix)\^-^dtdx 

for S small enough. Let E{x,t) = |P((9^/)(a:)p; because Pt is a smooth iden¬ 
tity with a convolution kernel it is elementary to show that sup|,j,_j,|<j|L"(j/,t)| < 
CM{d^f){x). Let 

dy{x,t) = ls^t^-i^/s\ef{A^(;e-y){x)\'^^dtdx. 

By Lemma l9.21 and the preceding remarks, we have that 

1 r rRR) 

/ / \M<c+ c\\^im,s- 

R \R\JrJo 
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This establishes the desired bound on 


/• 


9.6. Terms with a„+i = 0 and /3„+i > 0. We conclude this section by bounding 
©f 1 for multiindices (3 with Pn+i > 0. Recall that 

0f/(a;)= f Aa^(y)/(y)d2/. 

a„+i=0 

By a well known argument of Fefferman and Stein [FS72] . using decay of the 
kernel of ©f (that is, the bound (13.31) 1. we find that if k is large enough then 

||0fl||c < C + supi^ / / |©f(l4Q)(a;)| ——. 

Q \Q\ Jo Jq t 

Let {Fp)a := A^p I 4 Q; then Fp is an array-valued function. More precisely, 

let q be the number of multiindices C, G N" of length m; alternatively, q is the 
number of multiindices ( G of length m with Cn-i-i = 0. We will think of 

as the vector space of arrays of numbers indexed by such multiindices. Then for 
each /3, Fp is a function in 1 —>■ C^). 

Now, observe that 

(9.17) A|| div„,j| : ^ c«) L2(R" C«) 

is a bounded operator, where V™ and div„j y are as defined in Section [51 and so we 

have a Hodge decomposition of L^(]R"' 1 —>■ C^). Specifically, if F' G L^(]R" 1 —>■ C^), 
then 

F = A-hA||V™$ 

for some H G L^(R" !->■ C^) and some $ G kF^(M"), with div,„ y H = 0 and with 
II-H’IIl2(r"m.C'>) + II®IIw™>2(r"h^c) < C'll-P’lli2(RnM.co) • 

Applying the Hodge decomposition to Fp, we see that 

©fl4Q(x)= Y. [ dr^''~^d^E\x,t,y,0)[Hp + A\\VY^p) dy. 

a„+i=0 “ 


But because div,^,! y Hp = 0, we have that 


©fl 4 Q(x)= Y f dr^’^-^d:^E\x,t,y,0)(A\\Vf^p) dy. 
„„,1=0 V 


We may extend ^p to a function defined on by letting ^p{y,s) = ^piy). 

Observe that d^^p = 0 unless ^n-i-i = 0. Also, if a„+i = 0, then (AyV|p<i>^)a = 
{AV^^p)a. Thus, 

0fl4Q(a;) = Y X! / di;^^^~^d^E^{x,t,y,0)Aa(;{y)d^^p{y)dy. 

a„+i=0|C|=m 

Recall from formulas (14.31) and (12.311) that 

0f $/3)(a:) = - Y f d^~^^d^,E^{x,t,y,s)Aaciy)d^<3>p{y)dyds. 

|a| = |C|=m 
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Using the identity (12.2611 and integrating in s, we see that 

0f(Trm.| $/3)(a:) = t^f Ao^dy) dy. 

, , 1,1 JR" 


|a| = |C|=m 


Thus, 

0fl4Q(a;) 


= - E E^"/ dr^^-^d^,sEH^,t,y,0)A^^iy)d<^^iy)dy 

a„+i>0|C|=m 

+ 0 f(T>^,|d>^)(x). 


If an+i > 0, then a = 7 + e„+i for some multiindex 7 with I 7 I = m— 1. Conversely, 
if I 7 I = TO — 1, let 7 = 7 + Cn+i- We may write 

QtUQ{x)= E / ^r^''91sE^{x,t,y,0)A:y^{y)d‘^^fj{y)dy 

II 1 I /■ I 

|7|=m-l |C|=m 

+ 0f(Tr^,l$;3)(a;). 

By formula for 0f, we see that 

0 f l4Q(a;) = QtGp{x) + 0f (Tt^,| <i>f3){x) 


where (G/?)^, = Z^c 

Observe that Gp G L^(]R"). We may bound the term QfGp{x) as usual. 

Recall that 0f acts on the space WA"^ |{®"')5 the completion of {Tcm.l ‘fi ■ ‘fi & 
G^(]R"“''^)} under the norm. Consider the subspace lU, the completion of 

ip : ip & G“(]R"''''^), di^_^-yip{x, 0 ) = 0 for all a; G K" and all j > 1 } 

under the norm. We may let Ofdenote the restriction of Of* to the space W. 
Notice that <t>p = V|P$; 3 |r„, and so Of ^p){x) = 0f’"(Tt^j <^ 0 ){x). 

As we established in Sections 

||0f"l|k5<G + G||0fl||c.. 

and so we may control Of’^fTCm.l ‘h,s)( 2 ;)- This completes the argument that 0f 1 
satisfies a Carleson estimate. 


10. Test functions bg 

In this section we will choose test functions bg such that we may apply Theo¬ 
rem [ZTT] to bound 0f and 0f. (The remaining components of 0f were bounded 
in Section [5]) We will follow the example of |GdlHH] . which considers the case 
TO = 1. 

As in Section we will make the assumption 2m > n. Again, by Morrey’s 
inequality, this implies that functions locally in are locally Holder contin¬ 

uous. By Lemma |3.21 if 2m > n then solutions to elliptic equations are locally in 
L^(R” X { f|) for constants t, and thus are also locally Holder continuous. (See also 
AAA~*~lll Appendix B], in which a similar argument is made.) 

Fix some dyadic cube Q. Let yg be its midpoint. Let 

(10.1) Fs{x,t) = d'^~^E^{x,t,yg,s) 
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and let F± = F±i^i(^q-j for some small positive number k to be chosen later. By 
the bound (jQ.lOl) and the symmetry relation p.21|l . we may see that if n + 1 > 3 
then F±{x,t) G furthermore, by Theorem 12.201 we see that LF± = 0 in 

j^n+i i^y higher-order Green’s formula (I2.29L if t > 0 then 

drF-ix,t) = -drV^iTr+_,F.){x,t) + drS^(M+F_)ix,t) 
and by the corresponding formula in if t > 0 then 

0 = drV^iTr-_, F+)ix,t) + drS^iM^F+)ix,t). 

Adding and applying the definition (12.321) of P, we see that 

dTF_{x,t) = F+ - F_){x,t) 

+ dTs^{m\F_ + m-j,F+){x,t). 

Thus, by the definitions (14.1|) and (14.3|) of 0f and ©f*, 

(10.2) t>^d'^+^F_{x,t) 

= 0f F+ - F_){x) + 0f (M+ F_ + F+){x). 

Let 

(10.3) 6g = |Q|(T>„,|F+-T>„,|F_). 

Recall that is only defined as a linear functional on that 

is, as an operator acting on m — Ith-order traces of VL^-functions. Let bq be a 
representative of the operator -|- M(^F+); that is, bq is an array of 

functions that satisfies 

(10.4) (Trm-i = |Q|(Tfm-i Ma F- -1- T+)Rn 

= |0|(V’"^,AV™A_)^„+i + |Q|(V™(^, AV™F+)„„+. 

for all smooth, compactly supported functions (f. In Section [10.2l we will show that 
there is some such array of functions that in addition lies in L^(R."). 

Now, by formula (I2.26P and by definition of bq , bq and F _, 

&t{b^,b^q) = &?b^{x) + efbl^ix) = |Q|t'=ar+"F_(a;,t) 

= \Q\t'^drdT+'^-^E^{x, t, VQ, -Ki{Q)). 

An application of the bound (I9.10L with the roles of x and y reversed, reveals that 
the bound (I5.14|) is valid for this choice of bq = {bq, bq), albeit with constant Cq 
that depends on our choice of n. 

We thus need only show that this choice of bq satisfies the bounds ()5.15l) . ()5.16l) 
and (15.171) . with the distinguished component — bq in (15.161) the ej_-component 
of 

Remark 10.5. Although we will not make use of this fact, we observe that by the 
definition (I10.1|) of F±, the symmetry property (I2.21|) . and formula (I2.3ip for the 
double layer potential, we have that 

(5g,T>™_i(^)Rn = \Q\d!;:+^^V^\^m-iy^){yQ,-K£{Q)) 

- |Q|a™+YP^*(T>„_i ^){yq, KiiQ)). 
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Thus, bq may be viewed as the kernel of the double layer potential. If m = 1, 
then the classic jump relation Tr“ f — Tr"*^ / = / is well known. In the 

higher-order case, the analogous jump relation (see [Bar]) is 

Thus, if K is small enough and ip is smooth, then (bg, (^)Rn is approximately 

\Q\d^+i'^{yQ^^)- Thus, we expect bg to be approximately equal to e_L near Q, 
and so it is reasonable to expect the bounds (15.1611 and (|5.17|) to be valid for bg. 

Remark 10.6. Recall from formula (|2.25|) that if 2m > n -I- 1, precisely the case 
considered in this section, then the fundamental solution in the definition (110.111 
of Eg is only determined up to adding polynomials. However, note the presence of 
the vertical derivative in the definition of Eg] this vertical derivative suffices 

to remove the terms of the form fi^{x,t) (j/, s)'’ in formula (I2.25L leaving Eg{x,t) 
well-defined up to adding polynomials in x and t. The function Eg is a tool used 
to define bg and bg; notice from formulas (110.311 and (110.411 that these quantities 
depend only on the higher-order derivatives of Eg, and so the lower-order terms in 
formula (12.2511 do not affect our results. 


Remark 10.7. The conclusions of this section are also valid if n -I- 1 = 2; the 
analysis is somewhat more complicated because E± is no longer in 

By Morrey’s inequality, Lemma 13.21 and the bound (I2.23L we have that if 
n -I- 1 = 2 and i? is a cube of side-length \s — t\ then 


( 10 . 8 ) 



VZdT-'E^ix,t,y, 


s)pdy < 


|s-f| 


for any 5 > 0. In particular, 


[ \V^Eg{x,t)\^dx<-^. 

This is similar to the proof of the bounds (I3.3|l and (I9.10I1 . but we must use the 
bound (|2.23p instead of the bound (12.2211 in order to take m — 1 derivatives in the 
variable s rather than the variables {x,t). 

We may use the bound (19.1011 with the roles of x and y reversed to show that if 
t<s<aort>s>a, then 


f |V"^f,(x,t) - V^E^ix,t)\‘^dx < ^,*3 . 

Thus, if 0 < s < tr, then F±s — E±a- € 1T^(]R(^^^), and so we may apply the Green’s 
formula (12.2911 and the equivalent in to see that 


tkdrn+k 


F_s(x, t) - t) 


= ©t (Tcm.l Fs - Tr™,| F_g - Tr^j F^ + Tr^j F_^)(x, t) 
+ 0f (Mi F_g + M^Fg- M+ F,){x, t). 


Fix some t > 0 and let s = k£{Q). Observe that if we take the limit as cr —>■ 
oo, then the left-hand side approaches F-g{- ,t) in L^(K"). Furthermore, 

TCm.l-fio- —?> 0 in L^(]R"'). In Lemma 110.91 below, we will see that bg —>■ 0 as 
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the implied constant k 
^ 0 in L^{ 


as a 


oo; by definition of bg, this implies that F_cr + 
00 . Thus, by the bound (15.81) . we have that 


= 0f Fs -Tr^^iF_sK-,t) + 0f (M+ F_, + F,)( • , t) 

as L^(IR")-functions. Applying the Caccioppoli inequality, Lemma [XU and Mor- 
rey’s inequality, we see that this equality must be true pointwise as well. Thus, 
formula ( 110 . 21 ) is still valid if n + 1 = 2 and we may proceed as above. 

10.1. Bounds on bg. By the bounds (19.101) or (110.81) . we have that 


B" 






Thus, bg satisfies the bound (15.151) with constant Cq = Ck ”. 

We now show that bg satisfies the bound (15.171) . Following [GdlHHl Section 3], 
we fix a small positive constant u. Let (pg be supported on (1 + uj)Q with pg = 1 
on {1/2)Q. We may choose pg such that \V(j)Q{x)\ < 2f£{Q) for all x, and such 
that (pg > oj on Q. We then set d/ig = (pg dx. Observe that the conditions (j5.13l) 
are valid for Co = max( 2 , l/w). 

Then by definition of bg, if we let Fg = F+ — F_, then 


/tg 

JQ 


d^Q / Tr^ I Fg < 


[ '^m,\Fg(pg- ( Trrn,\Fg(pg. 
Jr" 9r"\q 


Recall that each component of Tr^_| Fg may be written as d^Fg{x, 0) for some (3 
with Pn+i < |,5| = TO. In particular, /3 = + 7 for some 1 < j < n and some 

multiindex 7 . Integrating by parts, we see that 


T^m,\ Fg (pg 


< 


|Trm-i Fg\ |V(/)q|. 


Recalling the regions on which 

c r 


and V^Q are supported, we see that 


£ 




dfJ-g 


< 




|V™-'Fg| + 


£ 


(l+u;)Q\Q 


|V"^Fq|. 


'(l+c^)Q\(l/2)Q 

Now, observe that LFg = 0 away from {yg, ±k£{Q)). Thus, we may apply Holder’s 
inequality and Lemma 13.21 to see that 


/ 

Jo 


bgdyg 


< 


C|Q|i/^ / r MQ) 

V2Q\(1/4)Q/-I!(g) 


^Fg\‘- 


1/2 


C 




riiQ)/2 \ 1/2 


\J{3/2)g\{i/2)gJ-e{g)/2 
We may use the Caccioppoli inequality (Lemma |33J to control the second term by 
the first term. For ease of notation let S = {2Q\{1/4)Q) x iiQ))- Recalling 

the definition of Fg, we see that 


jQ 


bQ 


dyg 


< 


C\Q\ 


rKe(Q) 


'-K.e(Q) 


Vyp^dTE^{x,t,yg,s)ds 


dx dt. 


2 
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Applying Holder’s inequality again, we see that 

/ b^d^XQ ' < \y':^^d:^E^{x,t,yQ,s)\^dsdxdt 

Using formula (12.261) . (12.211) and the bound (I9.10L we see that 

2 


{• 

JQ 


bg dfJ,Q 


< Ck^ 


Thus, if we choose k small enough, then 6 q satisfies the bound (15.171) . Notice that 
K may be chosen depending only on the constant Ci in the bound (I5.12L that is, 
on the numbers determined in Section 121 It is acceptable for the numbers Co in the 
bounds (15.131) . (15.141) and (15.151) to grow as k —>■ 0. In particular, recall that bg 
satisfies the bound (15.151) with a constant C'o(k) = this growth is acceptable. 

Thus, bg satisfies all the conditions of Theorem 15.Ill 


10.2. Bounds on bg. To conclude the proof of Theorem ll.61 at least in the case 
2m > n, we need only show that bg satisfies the bounds (15.151) . (15.161) and (15.171) . 

The most involved argument of this section will be the proof of the following 
lemma. 


Lemma 10.9. Suppose that 2m > n. If g G lid4^_]^(R"), then 

(10-10) l(9)bQ)R"l < 

Ki ' 

Furthermore, if g — 0 in (l/AfQ and n < 1/16, then we have a better estimate: 
(10.11) l(9,b|)Rn| < C K^/\Q\\\g\\L^-Rr,). 

The bound (jlO.lOl) is valid even for k large; recall that this bound was used in 
Remark 110.71 to show that the bound (15.141) is valid even in dimension n + 1 = 2. 

Notice that this implies that bg is a bounded linear functional on WA‘^_-^(RF); 
if m > 2 then this is a proper subspace of L^(R”). Thus, bg lies in a quotient 
space of L^(]R"). Once this lemma is proven we may extend bg to a bounded linear 
functional on L^(R”) (establishing the bound ()5.15|) b for example by orthogonal 
projection; we will need to select our projection carefully to ensure that bg, after 
projection, satisfies the bound (15.171) . 


Proof of Lemma \10.9[ It suffices to prove this lemma for all g such that g = 
Tr^-i rj for some smooth, compactly supported function g. Recall that 

(S) ^q)R" = IQKSj F_)Rn + \Q\{g, T’+)Rn 

= |Q|(V™G, + |Q|(V™G, AV”^F+)jj„+i 

for any extension G of g. 

We will need to construct our extension G of g carefully. Let H be the extension 
of g given by Lemma [SHI Recall that H satisfies the estimates (13.61) and (13.81) . and 
that if g = 0 in {\/A)Q then V^~^H{x,t) = 0 whenever \t\ < dist(a;,]R" \ {1/4)Q). 
Let p be smooth, supported in B{Q, 1/2) and integrate to 1. Suppose further that 
the higher moments are zero, that is, pix) dx = 0 for all 1 < |CI < w. Let 

G{x,t) = f H{z,t)dz = f ip{z)H{x-zt,t)dz. 
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To study the derivatives of G, observe that for some constants 


d'^G{x,t)= zt,t)dz 

iei=ici.cii>cii 


where C|| denotes the horizontal part of C, that is, Cy = (Ci,..., C„). Let J(^^^{z) = 


( 10 . 12 ) 





X — z 


d^H{z,t) dz. 


By our moment condition on (/?, we have that 



if C = ?, 


/ 

JR' 


t)dz = 0 


otherwise. 


Furthermore, is a smooth cutoff function, and so Trm-i G = g. So we need 
only bound 


|Q|(V™G, + |Q|(V™G, AV™F+)k„+i. 

We will need some special arguments to establish the bound (110.111) . Arguing as 
in the proof of Lemma [3.41 we see that if g = 0 in (l/4)(5 then V^~^G{x,t) = 0 
whenever 2|t| < dist(a:,M"' \ (l/4)(5). In particular, if k < 1/8 then V^^^G = 0 
near (j/g, ±kI{Q)). (We require k < 1/16 so that V™“^G = 0 everywhere within a 
fixed radius of (yg, ±kI{Q)).) Observe that LF+ = 0 away from these points, and 
so 


0 = (V™G, AV™F+)r„+i = (V™G, AV’"F+)„„+i + (V™G, AV™F+)g„+i. 

Thus, if g = 0 in (l/4)(5, we need only bound 

|Q|(V™G,AV™(A_-F+))r„+i. 

We now introduce some notation. Let G = {{x,t) ■ t > (1/2) dist(x, R" \ 
(l/4)(5)}, so if g = 0 in (l/4)(3 then supp V™“^G D C G- If |q;| = m, 
let 

wl^{x,t) = {A{x)V^Fs{x,t))a= X! Aap{x)d^Fs{x,t) 

|/3|=m 

with w'^ = Wa. Let Wa = — w~. To establish the bound (IIP.Ill) , we need 

only bound 



To establish the bound (|10.10p . we need only bound the quantity 

f d°‘G{x,t) w~{x,t) dx dt + ( d°‘G{xG)w'^{x,t)dxdt. 

I I Jr”+^ Jr"+i 

\a.\—m + - 

The second integral is similar to the first integral; thus, we will present the argument 
only for the first integral. In other words, we will work only in R"^^, not R"“''^, 
whether our goal is to establish the bound (IIP.101) or (110.111) . 
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We will need to bound w®; it will also help to bound vertical derivatives of w®. 
Let j > 0 be an integer. Observe that by formula (j2.26|) and the definition of Fs, 

dx = f dx. 

Jr^ 

By the bounds (19.1011 or pO.Sp . if 2m > n and j > 0 then 
f '^iw^^{x,t)f dx < C\t - s\-^-^F 


I 

Jr^' 


/R" 


Thus, we have the bounds 
(10.13) 


sup 

t>0 . 


' dx < 


(10.14) 

(10.15) 


C 

k”|Q| ’ 

c 


^ [ \dtw-{x,t)ftdtdx^ , 

R" Jo ^ lOI 

f \ 1/2 (J 


Now, by formula (12.2611 . d{w^{x,t) = {—iydiw^{x,t). Furthermore, if t > 0 
|s| < Ki{Q), then by the bound (19.1011 


/ 

Jr-^' 


lg{x,t) \dsW^{x,t)y dx < C{£(Q) +t) 


-n-2 


Thus, recalling that Wa{x,t) = w'^{x,t) — {x,t), we have that 


(10.16) 

(10.17) 


sup / lg{x, t) \Wa{x,t)y dx < 
t Jr" 

[ [ lgix,t)\dtWaix,t)\'^ tdtdx < 

7r" Jo IOI 

(S' \ 1/2 


IQI ’ 


(10.18) 

Recall that we wish to bound 

/ d°‘G{x,t)wa{x,t)dxdt or / d°‘G{x,t) w~{x,t) dx dt. 

We will essentially proceed by integrating by parts to move one derivative from 
G to Wa', we will need separate arguments in the case where we integrate by parts 
in t (possible only if a„+i > 0) and in the case where we integrate by parts in a 
horizontal variable Xj (possible only if a„+i < m). 

First, if an+i > 0, then 0 = 7 + e„+i for some multiindex 7 with I7I = m — 1. 
So 

I d°'G{x,t) w~{x,t) dx dt = I dtd"'G{x,t)w~{x,t) dx dt. 

Jr^+i iR”+^ 

Integrating by parts in t, we see that 

[ d°‘G{x,t)w~{x,t)dxdt = — lim f d''^G{x,t) w~{x,t) dx 

Jr'I + I t-S-0+ Jrti 
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Recall that H satisfies the uniform L? bound (18.6^ : by formula (I10.12|l . the same is 
true of G. We may control the first term using the estimate (110.1811 and the second 
term using the estimate (110.151) . This yields the bound 

d°‘G{x,t)w~{x,t)dxdt 


< 




l9l|L2(Rn). 


Similarly, 


/ d°‘G{x,t)wa{x,t)dxdt = — lim / d'^G{x,t) Wa{x,t) dx 

t^0+ 

— I d'^G{x,t) dtWa{x,t) dx dt. 

Jr"+i 

If g = 0 in (1/4)Q then we may integrate over Q rather than K"^^. By the bounds 
(110.181) and (110.161) on Wa and the uniform estimate on we have that 


/, 


d°‘G{x, t) Wa {x, t) dx dt < 


Gk 

7W\ 


l9llL2(Rn). 


Now, we turn to the case where a„+i < |a| = m. We still integrate by parts 
in t. We see that, if Wa = Wa or Wa = Wa, then 


/, 


d^G{x, t) Wa{x, t) dx dt 


-[ 

jR’i 


tdt{d°'G{x,t) Waix,t)^ dxdt 


td°‘dtG{x,t) Wa{xj t) dx dt — 


td°‘G{x,t) dtWa{xj t) dx dt. 


Recall that G as well as H satishes the estimate (|3.8ll , and so by the bounds (|10.14p 
and (|10.17l) . 

td°‘G{x,t)dtWa{x,t)dxdt ' ^ 


f 

7r^+i 


< 




l 9 'llL 2 (Rn), 


L 


t d°‘G{x, t) dtWa (x, t) dx dt 
We are left with the term 


GK II , I, 

< -^=l| 9 'llL 2 (Rn). 




/. 


td°'dtG{x,t)wa{x,t)dxdt, ci„+i < |a| = in. 


We have square-function estimates on dn+iWa rather than Wa] thus, we write 

f td°‘dtG{x,t)wa{x,t)dxdt = f td°‘dtG{x,t) f drWa{x,r) dr dx dt. 
Jr^+^ Jr'^+^ Jt 

Observe that if r > f > 0 and g = 0 in (l/4)(5, then lg{x,r) d°‘dtG{x,t) = 
d°‘dtG{x,t); this is true because, if lc;(a;,r) ^ 1, then V™G(a:,f) = 0. Let 
Va{x,r) = drWa{x,r) or lg{x,r) drWaix,r), depending on whether we seek to 
establish the bound (jlO.lOll or (IIO.IIL Then 

f td°dtG{x,t)waix,t)dxdt = f td°dtG{x,t) f Va{x,r) dr dx dt. 

Jr^+i Jr”+'- Jt 
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If an +1 < \a\ = m, let j = ja be any integer such that j < n and aj > 0, and let 
C = Ca = o: — Cj + e„+i. (For the sake of definiteness we may let ja be the smallest 
such integer.) Then 

f td°‘dtG{x,t)'Wa{x,t)dxdt = f tdj^d^°‘G{x,t) f Va{x,r)drdxdt. 

Jt 

For each pair of multiindices a and /3 with |a| = |/3| = m, define the linear operator 
Ta,i 3 by the relation 

Ta,l3F{x, 0 = ^ (“T~) 

where Ja,p is as defined above in the discussion of V^G and Then 

/ td°‘dtG{x,t)wa{x,t)dxdt 

= ^ H{x,t) f Va{x,r) dr dx dt. 

lai Jr”+^ Jt 

1 ^ 1=771 4 - 

We may rearrange the terms of the integral to see that 
f tdj^T(^^^ 0 d^H{x,t) f Va{x,r) dr dx dt 

’ Jt 

^ f f f f dxdr d^H{z,t) dzdt. 

Jo Jt Jr^ ^ \ ^ / 

We will use the Christ-Journe T1 theorem fTheorem 15.II above 1 to bound 
Wa,r{z,t) = J ^{dj^Jc^,l3)(^^^-j^'^Va{x,r) dx. 

Let ij}t(z,x) = (^^)- We then have that 

\Mz,x)\<^, \VxMz,x)\<-^ 

and both terms are zero if |a; — z| > t/2. Finally, observe that / ’tpt{z,x)dx = 0 
and so 0tl(a;) = 0; thus, we have the estimate 

pOO p 7 7i 

/ / \WaA^,t)\^ ^ <C\\Va{- ,r)\\l.(^„y 

Jo JR" ^ 


Thus, 


[ F[{x,t) f Va{x,r) dr dx dt 

jRj+^ ’ Jt 

poo ps p 

= / / Wa,riz,t)d^H{z,t)dzdtdr 

Jo Jo ^R" 

By the bound (j3.8|) and the above bound on Wa,r, 

pOO ps p poo 

/ / / Wa,r{z,i)d^H{z,t)dzdtdr < / ||ua( •, r)||i 2 (Rn) ||g||i 2 (Rn) dr 

Jo Jo Jr'^ Jo 
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and by the bounds (I10.15|) and (I10.18I1 . 



t d°‘dtG{x, t) {x, t) dx dt 


< 


C 




l9'llL2(Rn), 



t d°‘dtG{x, t) Wa {x, t) dx dt 


< 


Gk 

7W\ 


ll9lU2(Rn) 


as desired. 


□ 


We now have that bg is a bounded linear operator on a subspace of 

L^(K"). We extend bg using a similar projection as in Section ISTTl the difference in 
this case is that we use only two projection operators rather than countably many. 
Let Wn and Wf he the closure in of, respectively, 

Wn = {1(1/2)Q Ttm-l + (1 ” 1 (i/2)q)/ ^ V ^ G ^, f G L^(IR"')}, 

Wf = {1(1/4)q/ + (1 - 1(1/4)q) V^-.<pGG^,fG L2(R")}. 

Let On and Of denote orthogonal projection from onto the subspaces W„ 

and Wf, observe that Onf = f outside of {1/2)Q and that Off = / inside (l/djQ. 
Furthermore, 0„(Trm-i f) = 0/(Trm-i f) = T^m-i f for any nice function ip. 

Let T] be smooth, supported in {1/2)Q x {—^{Q)/A,^{Q)/A) and identically equal 
to 1 in (l/4)(5 X (—^(Q)/8,t!((5)/8), with the usual bounds on the derivatives of ry. 

Define 7r„ : Wn '->■ WAf^_^{M.'^) and ttj : Wf i-A as follows. Suppose 

that / = Trm_i in {1/2)Q or K" \ {1/A)Q for some smooth function tp. We 
may renormalize f so that /^•J^y2)Q\(l/4)Q Tr = 0 for all |C| < m — 1. Let 

TTnf = Trm-i(?7</5) and let tt// = Trm-i((l — v)^)- Observe that 7r„ and tt/ are 
well-defined, that 7r„/ = 0 outside {1/2)Q and that Tiff = 0 in {1/4)Q, and that by 
the Poincare inequality 7r„ : Wn >->■ and nf : Wf lW4^_j(R") are 

bounded operators. Finally, notice that 7r„(Trm-i f) -f 7r/(Trm_i (p) = Trm-i f 
for any smooth, compactly supported (p. 

We let bg satisfy 

(/) ^q)r" = {T^nOnf + TT/O/Z, 6§)Rn 

where the right-hand side is given by formula (110.411 . Notice that if / = Thm-i <p, 
then Onf = Off = f. Thus, this definition is consistent with formula (I10.4L 
By the bound (110.101) and boundedness of 7r„ and tt/, we see that 

I(/.^q)r”I < —p7^\\f\\L^{RfVW\ 

and so the bound (j5.15|) is established with Oq = Ck“". We are left with the 
bounds (jh.lOp and ()5.17p . 

Observe that by the bound (IIP.lip . 

|(7r/0//,f)§)Rn| < C'||/||L2(Rn)KV^. 

Furthermore, if / = 0 in (1/2)Q then iTnOnf = 0; thus, we have that 
(10.19) II^qIIl2(R"\(i/2 )Q) < Ck\/\Q\. 
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Fix some 7 with I 7 I = to — 1, and let bg = ( 6^)7 for some I 7 I = to — 1. Then 


bq = . We seek to show that 


Re 


f \ [ bl{x)(j)Q{x)dx > a or ^ f bl{x)(j)Q{x) dx 

Jq (PQ JQ Jq Jq 


< rja 


for some constant a independent of Q and some ij depending on Ci. Notice that 
lo ~ some constant c depending on (j)Q, with 1/2 < c < (1 + w)". 

Let 

= ^^ {x-yQ,ty(j)Q[x) p{t) 

where p{t) = 1 for |t| < ^{Q) and p{t) = 0 for |t| > 2(.{Q). 

Notice that if x G {1/2)Q and |t| < f’(Q), then V"''“^$q(x, t) = = (j)Q{x)e^. 

In particular, feg/iQ = Thm-i ■ f>Q in {1/2)Q. Furthermore, bgcjig =0 and 
Tr^-i ‘&Q ■ bg = 0 outside of (1 + uj)Q. Thus, by the bound (jl0.19L 


f bl{x)(j>Q{x)dx-[ Tr™_i$^-6Q < 7 ^ / 

\Q\ Jq ^ \Q\ Jr" ^ \Q\ J(i 


bgl < Ck. 


(l+u:)Q\(l/2)Q 

Thus, to establish the bounds (15.161) and (15.171) . it suffices to bound the quantity 


-r^ f T 

\Q\ Jr" 


Tr™_i 


from above or from below. 

Applying the definition of 6 q, we see that 

fr„_i • f>Q = v”^$^{x, t) ■ A{x)V^F_ (x, t) dx dt 

+ / V™$/,(x,t)-A(x)V"'F+(x,t)dxdt. 

Jr7+^ 

Now, observe that = 0 in {1/2)Q x {—J{Q), ^{Q)). Applying the definition 

of F± and the bounds (19.101) or (110.81) . we see that 


[ V™$X(x,t) • A(x)(V™F+(x,t) - V'”F_(x,t)) dxdt 

/r"+' 


< Ck 


and so we may consider 

f • AV^F. 

iR"+' ^ 

Now, recall that 


• AV'^F. = 


/ 

./R" + l 


• AV'^F.. 


f '^^<^>l{x,t) ■ AV^F_{x,t)dxdt 
/r^+i 

= [ V™$^(x,t)-A(x)V™ar-'if^(x,t,2/g,-«^(Q))dxdt. 

7r**+i 
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Applying the symmetry property (I2.21L we see that 


Rn + l 


{x, t) ■ (x, t) dx dt 


f dT ^V^^E^“‘{yQ,—K£{Q),x,t)-A*{x)\7"^^Q{x,t)dxdt. 

Jr"+i ’ ^ 


By formula (I2.24I1 . 


Rn + l 


t) ■ AV™F_(x, t) dx dt = (A*V™$^)( 2 /q, -k£(Q)). 


Recall (formula (|2.17|) 1 that 11^ (A*V'"<i>Q) = <i>Q. Thus 

V'"$^(x,<) • AW^F-{x,t)dxdt = ar-^$^( 2 /Q,-K€(Q)). 


The right-hand side is equal to one if 7 = 7_l and is zero otherwise, and so the 
bounds (15.161) and (15.171) are established. 


11. Reduction to operators of higher order 

We have now shown that 0^ and 0f satisfy the bounds (15.3p and (j5.4l) . We have 
established that whenever 2m > n, the condition (15.121) is valid, and there exist 
functions bg such that the conditions (15.141) . (I5.15|l . (15.161) and (I5.17|) are valid. 

Thus, if 2m > n, then by Theorem 15.111 and 0f satisfy the bound (j5.18|) : 
this implies that the bounds dni and m are valid. 

We now must establish these bounds for operators of order 2m < n. We use a 
fairly standard technique in the theory of higher-order differential equations; see 
[AHMTOH Section 2.2] and |Barl41 Section 5.4]. 

Fix some operator L of order 2m < n, and choose some number M such that 
2m -I- 4M > n. Now, there are constants ac^ such that 

A“= ^ 

|C|=M 

In fact, Oij = m!/^!, and so we have that a^- > 1 for all jCj = M. 

Define the differential operator L = A^LA^\ that is, (v?, Lij}) = {A^tp, LA^ip) 
for all nice test functions (p and !/>. We remark that L is associated to coefficients A 
that satisfy 

(11.1) A^^(^x) = 'y ] (x) = y ] —2C)(e—2^) {^)’ 

a+2(=S |C| = M, 2C<5 

l3+2i=e \i\=M, 2{<e 


Observe that A is t-independent and satishes the bounds (12.41) and (12.51) . It was 
shown in the proof of |Barl41 Theorem 62] that 

E^{x,t,y,s)= a<:;a^dl‘^tdl^^E^{x,t,y,s). 

ICI = ICI=M 
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Now, by formula (I2.35L if |a| = m then 
d°‘S^g{x,t)= [ d^^tdl,E^{x,t,y,0)gj{y)dy 

|7|=m-l|C|=M|4|=M 

Let ge{x) = Y.-y+ 2 i=e‘^^ 9 j{x). Notice that \g{x)\ < C\g{x)\. Then 
(11.2) d°‘S^g{x,t)= Y ac^d°‘^'^‘'S^g{x,t). 

|C|=M 

Thus, because the bound (ini is valid for operators L of order 2m + 4M for M 
large enough, we have that 

p /»00 p pOO __ 

/ / \'^^dtS^g{x,t)ftdtdx<C / \V"^+'^^dtS^g{x,t)ftdtdx 

Jr-^ Jo Jr^ Jo 

^ C'llfl'llL2(][{n) < C'||9|lL2(J{n) 

and so the bound uni) is valid even for operators of order 2 m < n. 

The argument for 2?^ is somewhat more involved. In this case we will use 
Theorem l5.2l observe that ©f satisfies the bounds (15.31) and (15.4L and so we need 
only establish the bound (15.5p . that is, to bound ©fe^ for multiindices /3 with 
\I3\ = m. 

Recall from formula (19.11) that 

Qfefi{x) = - X ^E^{x,t,y,s) Aafi{y)dsdy 

I I Jk^+i 

|Q;|=m — 

and so 

epix) = - Y X sE^{x,t,y,s)Bsp{y)dsdy 

|C| = M |<5|=m+2M 

where 

Bspiy) = Y «C^a/3(y)= X A<5-2«)/3(2/)- 

a+2{=5 |5|=M. 2J<<5 

Recall our formula (111.11) for the coefficients A ot L. We then have that 

Aseiy) = X ««^<5(e-2o(2/)- 

|{|=m. 25 <e 

We claim that there exist numbers such that 

Es^ — 'y ^ ^/5e A§^. 

\€\—m-\-2M 

To see this, define the operator iL by = X]|j|=m 2 ^<e ^e- 2 {- It suffices to 

show that we may recover F from 4i(F). Begin with indices a such that aj < I 
for all but one value jo of j. Let ^ = Mej^ and let e = a + 2^. Then ^ is the only 
multiindex with |^| = M and with 2^ < e. Thus 4'(F')e = Fa = Fa', in other 
words, we can recover Fa from 'L(F'). Next, consider indices a with aj^ arbitrary, 
2 < cijj <3, and aj < 1 for all other values of j. Let ^ be as before and let 
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e = a + 2^. Then 'P(-F)e = Fa + a( Fe- 2 C, where C = eji + [M — l)ejg. Since 
Fs- 2 q is known we may recover Fa- Continuing in this fashion, we may recover Fa 
from il'(j') for all multiindices a. 

Thus 


0fe 


= ^ bpet^ f^^dr^^dl‘^tdy^^E^{x,t,y,s)Aseiy)dsdy. 

\C\=M \S\=m+2M 

\e\—m-\-2M 

By formula (12.311) for V and (12.321) for V, this equals 

efe0(x)= dl‘^tD^ee{x,t)- 

|(^|=M |e|=m+2M 

Recall from formula (|4.3I) that j{x) = Define 

0f/(a:)=t'='ar+"^+'='5^/(a:,t) 

for some k' to be chosen momentarily. 

Because m + 2M > n, if k' is large enough then we have that satisfies the 
estimates (15.181) and (15.3L and so by Lemma lOl 

r-dQ) 


1 


sup ■ 
Q 


fc'am+2M+fc' j-iA 


\t^d: 


Q JO 


D^es{x,t)\' 


, dt dx 
t 


< C. 


Fix some cube Q and observe that 
r^Q) 


/ / \Q?ep{x)\ 

JO JQ 


9 dx dt ^ V—^ 

<c ^ 


t 


\€\—m-\-2M 


e(Q) 


0 JQ 


Jo 


Applying the Caccioppoli inequality in Whitney boxes, we see that 


FiQ) 

Jo 


|0t e/3(T)|' 


, dx dt 
t 


E 

|£|=m+2M 


2dQ) r 
J 2 Q 


dx dt 


, dx dt 


0 JQ 

If we let k = 2M F k', we see that 

rt{Q) r fjf 

/ / \Qfep{x)\^^<C\Q\ 

JO JQ * 

and so 0^* satisfies the bound (15.51) . Thus, by Theorem l5.2l we have that 0^ satisfies 
the bound (15.61) . Thus, by Lemma I4l5] we have that satisfies the bound (II.8p . 
as desired. 


References 

[AAll] Pascal Auscher and Andreas Axelsson, Weighted maximal regularity estimates and 

solvability of non-smooth elliptic systems J, Invent. Math. 184 (2011), no. 1, 47-115. 
MR 2782252 El 

[AAA'^11] M. Angeles Alfonseca, Pascal Auscher, Andreas Axelsson, Steve Hofmann, and Seick 
Kim, Analyticity of layer potentials and solvability of boundary value problems 
for divergence form elliptic equations with complex L°° coefficients, Adv. Math. 226 
(2011), no. 5, 4533-4606. MR 2770458 [3lfil [Mlfikl 
[AAH08] Pascal Auscher, Andreas Axelsson, and Steve Hofmann, Functional calculus of Dime 
operators and complex perturbations of Neumann and Dirichlet problems, J. Funct. 

Anal. 255 (2008), no. 2, 374-448. MR 2419965 (2009h:35079) El El 














SQUARE FUNCTION ESTIMATES ON LAYER POTENTIALS 


59 


[AAMIO] 

[Agm57] 

[Agr07] 

[Agr09] 

[AHMTOl] 

[AM14] 

[AP98] 

[AQOO] 

[AR12] 

[AT98] 

[Bar] 

[Baris] 

[Barl4] 

[BM] 

[BMlSa] 

[BMlSb] 

[BM15] 

[CamSO] 

[CFK81] 

[CG83] 


Pascal Auscher, Andreas Axelsson, and Alan McIntosh, Solvability of elliptic sys¬ 
tems with square integrable boundary data, Ark. Mat. 48 (2010), no. 2, 253-287. 
MR 2672609 (2011h:35070) |1[5] 

Shmuel Agmon, Multiple layer potentials and the Dirichlet problem for higher order 
elliptic equations in the plane. I, Comm. Pure Appl. Math 10 (1957), 179-239. 
MR 0106323 (21 #5057) [7] 

M. S. Agranovich, On the theory of Dirichlet and Neumann problems for linear 
strongly elliptic systems with Lipschitz domains, Funktsional. Anal, i Prilozhen. 41 
(2007), no. 4, 1-21, 96, English translation: Funct. Anal. Appl. 41 (2007), no. 4, 
247-263. MR 2411602 120096:350701 [3l l6l lUl [Tsi 

_, Potential-type operators and conjugation problems for second-order strongly 

elliptic systems in domains with a Lipschitz boundary, Funktsional. Anal, i Prilozhen. 
43 (2009), no. 3, 3-25. MR 2583636 (20116:35362)0 

Pascal Auscher, Steve Hofmann, Alan McIntosh, and Philippe Tchamitchian, The 
Kato square root problem for higher order elliptic operators and systems on R", 
J. Evol. Equ. 1 (2001), no. 4, 361-385, Dedicated to the memory of Tosio Kato. 
MR 1877264 l2003a:35046') 151 155115511551 

Pascal Auscher and Mihalis Mourgoglou, Boundary layers, Rellich estimates and ex¬ 
trapolation of solvability for elliptic systems, Proc. Bond. Math. Soc. (3) 109 (2014), 
no. 2, 446-482. MR 3254931011 

Vilhelm Adolfsson and Jill Pipher, The inhomogeneous Dirichlet problem for 
in Lipschitz domains, J. Funct. Anal. 159 (1998), no. 1, 137-190. MR 1654182 
(99m:35048) [18] 

P. Auscher and M. Qafsaoui, Equivalence between regularity theorems and heat kernel 
estimates for higher order elliptic operators and systems under divergence form, J. 
Funct. Anal. 177 (2000), no. 2, 310-364. MR 1795955 (2001j:35057) [19] 

Pascal Auscher and Andreas Rosen, Weighted maximal regularity estimates and solv¬ 
ability of nonsmooth elliptic systems, II, Anal. PDE 5 (2012), no. 5, 983-1061. 
MR 3022848 0 

Pascal Auscher and Philippe Tchamitchian, Square root problem for divergence 
operators and related topics, Asterisque (1998), no. 249, viii+172. MR 1651262 
(2000c:47092) [36] 

Ariel Barton, Layer potentials for higher-order elliptic systems with rough coeffi¬ 
cients, in preparation. 1141 11811471 

_, Elliptic partial differential equations with almost-real coefficients, Mem. 

Amer. Math. Soc. 223 (2013), no. 1051, vi+108. MR 3086390 0 0 

_, Gradient estimates and the fundamental solution for higher-order elliptic 

systems with rough coefficients, ArXiv e-prints 120141. [5l 1141 [TSi 1161 [T9l 1561 
Ariel Barton and Svitlana Mayboroda, Layer potentials and boundary-value problems 
for second order elliptic operators with data in Besov spaces, Mem. Amer. Math. 
Soc., to appear. mm 

_, The Dirichlet problem for higher order equations in composition form, J. 

Funct. Anal. 265 (2013), no. 1, 49-107. MR 3049881 [S]H [6] 

Kevin Brewster and Marius Mitrea, Boundary value problems in weighted Sobolev 
spaces on Lipschitz manifolds, Mem. Differ. Equ. Math. Phys. 60 (2013), 15—55. 
MR 3288169 [3 [18] 

Ariel Barton and Svitlana Mayboroda, Boundary-value problems for higher-order 
elliptic equations in non-smooth domains, Harmonic Analysis and Operator Theory: 
Cora Sadosky Memorial Seminar in Analysis 1 (2015), to appear. (3 
S. Campanato, Sistemi ellittici in forma divergenza. Regolaritd aWinterno, Quader- 
ni. [Publications], Scuola Normale Superiore Pisa, Pisa, 1980. MR 668196 (831:35067) 
[19] 

Luis A. Caffarelli, Eugene B. Fabes, and Carlos E. Kenig, Completely singular 
elliptic-harmonic measures, Indiana Univ. Math. J. 30 (1981), no. 6, 917-924. 
MR 632860 (83a:35033) [3] 

Jonathan Cohen and John Gosselin, The Dirichlet problem for the biharmonic equa¬ 
tion in a domain in the plane, Indiana Univ. Math. J. 32 (1983), no. 5, 635-685. 
MR 711860 (85b:31004)[3 






60 

[CG85] 

[Chr90] 

[CJ87] 

[CMM82] 

[Dav95] 

[DJ84] 

[DJS85] 

[DK87] 

[DKV86] 

[DV90] 

[Eva98] 

[FJK84] 

[FMM98] 

[FS72] 

[GdlHH] 

[HKMP15a] 

[HKMP15b] 

[HMM13] 

[HMM15] 

[HoflO] 


ARIEL BARTON, STEVE HOFMANN, AND SVITLANA MAYBORODA 


_, Adjoint boundary value problems for the biharmonic equation on do¬ 
mains in the plane, Ark. Mat. 23 (1985), no. 2, 217-240. MR 827344 (88d:31006) [6] 

El [13] 

Michael Christ, A T{b) theorem with remarks on analytic capacity and the Cauchy 
integral, Colloq. Math. 60/61 (1990), no. 2, 601-628. MR 1096400 (92k:42020) [27] 
Michael Christ and Jean-Lin Journe, Polynomial growth estimates for multilinear 
singular integral operators, Acta Math. 159 (1987), no. 1-2, 51—80. MR 906525 
(89a:42024)[8][23] 

R. R. Coifman, A. McIntosh, and Y. Meyer, L’integrale de Cauchy definit un 
operateur borne sur pour les courbes lipschitziennes, Ann. of Math. (2) 116 
(1982), no. 2, 361-387. MR 672839 (84m:42027) [H 

E. B. Davies, Uniformly elliptic operators with measurable coefficients, J. Funct. 
Anal. 132 (1995), no. 1, 141-169. MR 1346221 (97a:47074) Eg] 

Guy David and Jean-Lin Journe, A boundedness criterion for generalized Calderon- 
Zygmund operators, Ann. of Math. (2) 120 (1984), no. 2, 371-397. MR 763911 
(85k:42041)[23] 

G. David, J.-L. Journe, and S. Semmes, Operateurs de Calderon-Zygmund, fonctions 
para-accretives et interpolation. Rev. Mat. Iberoamericana 1 (1985), no. 4, 1-56. 
MR 850408 (88f:47024) [27] 

Bjorn E. J. Dahlberg and Carlos E. Kenig, Hardy spaces and the Neumann problem 
in L'P for Laplace’s equation in Lipschitz domains, Ann. of Math. (2) 125 (1987), 
no. 3, 437-465. MR 890159 (88d:35044) [3] 

B. E. J. Dahlberg, C. E. Kenig, and G. C. Verchota, The Dirichlet problem for the 
biharmonic equation in a Lipschitz domain, Ann. Inst. Fourier (Grenoble) 36 (1986), 
no. 3, 109-135. MR 865663 (88a:35070) [6] 

Bjorn E. J. Dahlberg and Greg Verchota, Calerkin methods for the boundary integral 
equations of elliptic equations in nonsmooth domains. Harmonic analysis and partial 
differential equations (Boca Raton, EL, 1988), Contemp. Math., vol. 107, Amer. 
Math. Soc., Providence, RI, 1990, pp. 39-60. MR 1066469 (91i:35060) [1] 

Lawrence C. Evans, Partial differential equations. Graduate Studies in Mathemat¬ 
ics, vol. 19, American Mathematical Society, Providence, RI, 1998. MR 1625845 
(99e:35001)El][l9l[37] 

Eugene B. Fabes, David S. Jerison, and Carlos E. Kenig, Necessary and sufficient 
conditions for absolute continuity of elliptic-harmonic measure, Ann. of Math. (2) 
119 (1984), no. 1, 121-141. MR 736563 (85h:35069) [4] 

Eugene Fabes, Osvaldo Mendez, and Marius Mitrea, Boundary layers on Sobolev- 
Besov spaces and Poisson’s equation for the Laplacian in Lipschitz domains, J. 
Funct. Anal. 159 (1998), no. 2, 323-368. MR 1658089 (99j:35036) [3| 

C. Fefferman and E. M. Stein, spaces of several variables, Acta Math. 129 (1972), 
no. 3-4, 137-193. MR 0447953 (56 #6263) [44] 

Ana Grau de la Herran and Steve Hofmann, A local Tb theorem with vector-valued 
testing functions, preprint. [5] [S] 1241 1251126113011331135114511481 

Steve Hofmann, Carlos Kenig, Svitlana Mayboroda, and Jill Pipher, The regularity 
problem for second order elliptic operators with complex-valued bounded measurable 
coefficients. Math. Ann. 361 (2015), no. 3-4, 863-907. MR 3319551 [3] [4] 

_, Square function/non-tangential maximal function estimates and the Dirich¬ 
let problem for non-symmetric elliptic operators, J. Amer. Math. Soc. 28 (2015), 
no. 2, 483-529. MR 3300700 H [5] 

Steve Hofmann, Marius Mitrea, and Andrew Morris, The method of layer potentials 
in LP and endpoint spaces for elliptic operators with L°° coefficients, ArXiv e-prints 
(2013). [3] [5] 

Steve Hofmann, Svitlana Mayboroda, and Mihalis Mourgoglou, Layer potentials and 
boundary value problems for elliptic equations with complex L°° coefficients satisfy¬ 
ing the small Carleson measure norm condition, Adv. Math. 270 (2015), 480—564. 
MR 3286542 [3] [5] 

Steve Hofmann, Local T(b) theorems and applications in PDE, Harmonic analysis 
and partial differential equations, Contemp. Math., vol. 505, Amer. Math. Soc., 
Providence, RI, 2010, pp. 29-52. MR 2664559 f2011e:42024) [771 




SQUARE FUNCTION ESTIMATES ON LAYER POTENTIALS 


61 


[JK81] 

[KKPTOO] 

[KP93] 

[KR09] 

[KSll] 

[MitOS] 

[MM85] 

[MMll] 

[MM13a] 

[MM13b] 

[MMSIO] 

[MMWll] 

[Nad63] 

[PV92] 

[PV95a] 

[PV95b] 

[Rosl3] 

[Rul07] 

[Sem90] 

[She06a] 

[She06b] 


David S. Jerison and Carlos E. Kenig, The Dirichlet problem in nonsmooth domains, 
Ann, of Math. (2) 113 (1981), no. 2, 367-382. MR 607897 (84j:35076)|l 
C. Kenig, H. Koch, J. Pipher, and T. Toro, A new approach to absolute continuity 
of elliptic measure, with applications to non-symmetric equations, Adv. Math. 153 
(2000), no. 2, 231-298. MR 1770930 (2002f:35071) H 

Carlos E. Kenig and Jill Pipher, The Neumann problem for elliptic equations with 
nonsmooth coefficients. Invent. Math. 113 (1993), no. 3, 447—509. MR 1231834 
(95b:35046) H 

Carlos E. Kenig and David J. Rule, The regularity and Neumann problem for non- 
symmetric elliptic operators, Trans. Anier. Math. Soc. 361 (2009), no. 1, 125—160. 
MR 2439401 (2009k: 35050) |3] H] 

Joel Kilty and Zhongwei Shen, The LP regularity problem on Lipschitz domains, 
Trans. Amer. Math. Soc. 363 (2011), no. 3, 1241-1264. MR 2737264 (2012a:35072) 

[2] [18] 

Dorina Mitrea, A generalization of Dahlberg’s theorem concerning the regularity of 
harmonic Green potentials, Trans. Amer. Math. Soc. 360 (2008), no. 7, 3771—3793. 
MR 2386245 (2009b:35046) Ej 

Alan McIntosh and Yves Meyer, Algebres d’operateurs definis par des integrales sin- 
gulieres, C. R. Acad. Sci. Paris S&. I Math. 301 (1985), no. 8, 395-397. MR 808636 
(87b:47053) [27] 

Dorina Mitrea and Irina Mitrea, On the regularity of Green functions in Lipschitz do¬ 
mains, Comm. Partial Differential Equations 36 (2011), no. 2, 304—327. MR 2763343 

m 

Irina Mitrea and Marius Mitrea, Boundary value problems and integral operators for 
the bi-Laplacian in non-smooth domains, Atti Accad. Naz. Lincei Rend. Lincei Mat. 
Appl. 24 (2013), no. 3, 329-383. MR 3097019 El [7] 

_, Multi-layer potentials and boundary problems for higher-order elliptic sys¬ 
tems in Lipschitz domains, Lecture Notes in Mathematics, vol. 2063, Springer, Hei- 
delberg, 2013. MR 3013645 [2l l6l 171 [Til [Tsl 

V. Maz’ya, M. Mitrea, and T. Shaposhnikova, The Dirichlet problem in Lipschitz 
domains for higher order elliptic systems with rough coefficients, J. Anal. Math. 110 
(2010), 167-239. MR 2753293 ('2011m:350881 I^ITsl 

I. Mitrea, M. Mitrea, and M. Wright, Optimal estimates for the inhomogeneous 
problem for the bi-Laplacian in three-dimensional Lipschitz domains, J. Math. Sci. 
(N. Y.) 172 (2011), no. 1, 24-134, Problems in mathematical analysis. No. 51. 
MR 2839870 (2012h:35056) [TS] 

A. Nadai, Theory of Flow and Fracture of Solids, vol. II, McGraw-Hill, 1963. fT^ 

Jill Pipher and Gregory Verchota, The Dirichlet problem in L^ for the bihar¬ 
monic equation on Lipschitz domains, Amer. J. Math. 114 (1992), no. 5, 923—972. 
MR 1183527 (94g:35069) [3l[6] 

J. Pipher and G. C. Verchota, Maximum principles for the polyharmonic equation on 
Lipschitz domains, Potential Anal. 4 (1995), no. 6, 615—636. MR 1361380 (961:35021) 

[2] [18] 

Jill Pipher and Gregory C. Verchota, Dilation invariant estimates and the boundary 
Carding inequality for higher order elliptic operators, Ann. of Math. (2) 142 (1995), 
no. 1, 1-38. MR 1338674 (96g:35052) 0 [18] 

Andreas Rosen, Layer potentials beyond singular integral operators, Publ. Mat. 57 
(2013), no. 2, 429-454. MR 3114777 [3][5] 

David J. Rule, Non-symmetric elliptic operators on bounded Lipschitz domains in 
the plane. Electron. J. Differential Equations (2007), No. 144, 1-8. MR 2366037 
(2008m:35070)[3l[l[5] 

Stephen Semmes, Square function estimates and the T{b) theorem, Proc. Amer. 
Math. Soc. 110 (1990), no. 3, 721-726. MR 1028049 (91h:42018) [27] 

Zhongwei Shen, The L^ Dirichlet problem for elliptic systems on Lipschitz domains. 
Math. Res. Lett. 13 (2006), no. 1, 143-159. MR 2200052 (20071:35067) [2] [TS] 

_, Necessary and sufficient conditions for the solvability of the L^ Dirichlet 

problem on Lipschitz domains. Math. Ann. 336 (2006), no. 3, 697—725. MR 2249765 
f2008e:35059)f^rTHI 




62 

ARIEL BARTON, STEVE HOFMANN, AND SVITLANA MAYBORODA 

[She07a] 

. The. IT hoii.ndarii uahi.e problem.s on TAp.'=tr.hitz dnm.ains. Adv. Math. 216 
(2007), no. 1, 212-254. MR 2353255 (2009a:35064) [T3l 

[She07b] 

. A relationship betiaeen the Dirir.hlet and. regidn.rif.ii prohlem.s for elliptic, 
equations, Math. Res. Lett. 14 (2007), no. 2, 205-213. MR 2318619 (2008c:35043) [6] 

[Ste93] 

Elias M. Stein, Harmonic analysis: real-variable methods, orthogonality, and oscil¬ 
latory integrals, Princeton Mathematical Series, vol. 43, Princeton University Press, 
Princeton, NJ, 1993, With the assistance of Timothy S. Murphy, Monographs in 
Harmonic Analysis, III. MR 1232192 ('95c:42002) 1431 

[Swe09] 

Guido Sweers, A survey on boundary conditions for the biharmonic, Complex Var. 
Elliptic Equ. 54 (2009), no. 2, 79-93. MR 2499118 (2010d:35065) [12] 

[Tri78] 

Hans Triebel, Interpolation theory, function spaces, differential operators, North- 
Holland Mathematical Library, vol. 18, North-Holland Publishing Co., Amsterdam, 
1978. MR 503903 (80i:46032b) EH 

[Tri83] 

. Theory of fii.nr.tion .‘Spares. Monop-raphs in Mathematics, vnl. 78, Rirkhanser 

Verlag, Basel, 1983. MR 781540 (86j:46026) ED 

[Ver84] 

Gregory Verchota, Layer potentials and regularity for the Dirichlet problem for 
Laplace’s equation in Lipschitz domains, J. Funct. Anal. 59 (1984), no. 3, 572—611. 
MR 769382 (86e:35038) 0 

[Ver90] 

. The Dirichlet problem, for the poliiharm.onir equation in TA.psrh.i.t.z d.om.n.ins, 

Indiana Univ. Math. J. 39 (1990), no. 3, 671-702. MR 1078734 (91k:35073) ElllIIll 

[Ver96] 

Gregory C. Verchota, Potentials for the Dirichlet problem in Lipschitz domains, 
Potential theory—ICPT 94 (Kouty, 1994), de Gruyter, Berlin, 1996, pp. 167-187. 
MR 1404706 ('971:35041'! [2ll8l[T8l 

[Ver05] 

, The hi.h.n.rm.oni.r. Ne.um.n.nn problem, in Td.p.<^r.h.it.z d.om.n.in.<^. Acta Math. 194 
(2005), no. 2, 217-279. MR 2231342 (2007d:35058'l [6l 171 [TSl 

[Verio] 

. Ronnd.a.ry coe.rci.ve.ne..‘is n.nd. the. Ne.um.n.nn problem, for ff.h order l.ine.n.r pa.r- 
tial differential operators, Around the research of Vladimir Maz’ya. H, Int. Math. Ser. 
(N. Y.), vol. 12, Springer, New York, 2010, pp. 365-378. MR 2676183 (2011h:35065) 

m 

[ZanOO] 

Daniel Z. Zanger, The inhomogeneous Neumann problem in Lipschitz domains, 
Comm. Partial Differential Equations 25 (2000), no. 9-10, 1771—1808. MR 1778780 
(2001g:35056) El 


Ariel Barton, 202 Math Sciences Bldg., University of Missouri, Columbia, MO 65211 
E-mail address: bartonaeSmissouri. edu 

Steve Hofmann, 202 Math Sciences Bldg., University of Missouri, Columbia, MO 
65211 

E-mail address: hofmaiins@missouri.edu 

Svitlana Mayboroda, Department of Mathematics, University of Minnesota, Min¬ 
neapolis, Minnesota 55455 

E-mail address: svitlaiiaamath.umn.edu 








